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PREFACE 


This  final  report  on  "Strain  Induced  Waveguiding  in  Optical  Materials  for 
Integrated  Optics  Application",  ONR  Contract  Number  N00014-75-C-1097, 
covers  the  research  progress  made  and  results  achieved  during  the  period 
from  May  1,  1975  to  September  30,  1976.  This  work  was  carried  out  at 
Honeywell  Corporate  Research,  under  the  monitoring  direction  of  Lt.  David 
G.  Janssen  and  Mr.  J.  Trimble  of  the  Mathematical  and  Information  Science 
Division  of  the  Office  of  the  Naval  Research. 

Significant  research  results  achieved  under  this  investigation  are  summarized 
as  follows: 

• Techniques  for  preparing  laser  heating  strain- induced  (LHSI) 
optical  waveguides  were  developed,  and  the  properties  of  wave- 
guides made  with  optical  filter  glass  were  obtained. 

• A microprocessor-controlled  LHSI  optical  waveguide  circuit 
fabrication  facility  was  constructed,  and  simple  waveguide 
circuit  elements  were  obtained. 

• X-ray  and  optical  measurement  techniques  were  developed  for 
the  characterization  of  the  photoelastic  effect  and  the  effect  of 
optical  polishing  in  selected  integrated  optic  materials. 

Acknowledgement: 
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SECTION  I 
INTRODUCTION 


T echniques  of  guiding  optical  waves  in  planar  optical  materials  or  cylindrical 
optical  fibers  have  been  advancing  rapidly  in  recent  years.  The  major  objec- 
tive of  this  research  is  to  explore  the  micro-optic  fiber-optic  approach  for 
optical  communication  and  data  processing  applications,  with  the  ultimate 
goal  of  achieving  integrated  optical  devices  similar  to  the  integrated  elec- 
tronics. 

In  the  course  of  the  integrated  optics  research  it  was  observed  that  mechani- 
cal stress  and  strain  can  have  profound  effect  to  the  waveguide  properties  and 
device  performance.  Stress  and  strain  can  be  induced  by  lattice  mismatch 
in  epitaxial  growth  of  material,  by  thermal  expansion  or  contraction,  by 
inclusion  of  impurities  as  in  the  diffusion  operation,  and  by  grinding  and 
polishing  of  the  material  surface.  In  fact,  it  is  safe  to  assume  that  all 
material  surfaces  exhibit  some  degree  of  the  strain  effect.  These  residual 
or  externally  induced  strains  modify  the  waveguiding  properties  through  the 
photoelastic  or  possibly  other  effects.  For  example,  in  an  earlier  ONR- 
supported  research  effort,  we  have  demonstrated  that  optical  polishing  on 
PLZT  ceramic  can  cause  sufficient  refractive  index  increase  on  the  material 
surface  to  result  in  waveguiding^. 

Besides  this  interesting  example,  it  is  expected  that  similar  effect  can  be 
found  in  other  integrated  optic  material  where  the  surface  is  prepared  by 
optical  polishing.  The  technique  to  examine  the  stress  and  strain  on  an 
optically  polished  surface,  and  the  correlation  of  polishing  induced  strain 
to  the  waveguiding  properties  on  integrated  optic  media  are  important  tasks 
of  integrated  optics  research. 
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Alternatively,  an  attractive  technique  for  the  fabrication  of  integrated  optic  » 

multimode  channel  waveguide  is  to  use  laser  heating  to  introduce  localized 
strain  in  an  otherwise  uniform  medium.  This  approach  is  adaptable  for  com- 
puter control  to  generate  the  desired  multimode  channel  waveguide  patterns. 

Realizing  the  importance  of  the  effect  of  strain  on  integrated  optic  applica- 
tions, the  present  research  was  initiated  by  the  Mathematical  and  Information 
Science  Division  of  the  Office  of  Naval  Research,  and  a contract  was  awarded 
to  the  Honeywell  Corporate  Research  Center  to  perform  this  research  study. 

The  first  phase  of  the  present  contract  began  from  May  1 and  ended  Decem- 
ber 31,  1975.  A follow-on  contract  was  granted  to  cover  the  period  ending 
September  30,  1976,  for  the  initial  work  of  the  second  phase  research. 

The  objectives  of  this  research  are: 

• To  determine  the  usefulness  of  the  strain- induced  effect  for  inte- 
grated-optics applications.  The  strain  induced  refractive  index 
change  on  suitable  integrated  optics  materials  will  be  of  primary 
concern. 

• To  obtain  basic  understanding  and  perform  measurements  of  the 
effects  of  strain  on  the  variou.-  physical  phenomena  in  optical 
materials  useful  for  integrated  optics  applications. 

• To  develop  device  concepts  using  selected  materials  and  methods. 

The  use  of  laser  beam  heating  for  strain- induced  integrated  optic 
circuit  will  be  emphasized. 

These  objectives  have  guided  our  research  during  this  period  of  investiga- 
tion. Specific  tasks  performed  during  this  research  effort  are: 
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m To  develop  techniques  of  inducing  controlled  residual  strain  for 
optical  waveguiding. 

• To  measure  the  effect  of  strain  on  selected  materials  and  to 
correlate  them  to  the  properties  relevant  to  optical  waveguiding. 

• To  produce  strain-induced  planar  waveguide  or  channel  waveguide 
and  obtain  the  waveguide  properties. 

• To  fabricate  simple  channel  waveguide  circuits  by  laser  heating. 

• To  construct  a microprocessor  controlled  facility  for  automated 
fabrication  of  laser- induced  channel  waveguides. 

In  this  final  report,  we  present  the  progress  made  and  results  obtained  in  all 
the  above  areas.  The  significant  progress  and  results  are: 

• Techniques  for  preparing  laser  heating  strain- induced  (LHSI) 
waveguides  were  obtained  and  the  waveguide  properties  measured. 

• Techniques  for  automated  fabrication  of  LHSI  waveguide  circuits 
were  developed  and  a microprocessor  controlled  fabrication 
facility  constructed. 


• X-ray  and  optical  techniques  were  developed  for  the  characteriza- 
tion of  residual  or  induced  strain  in  integrated  optic  materials. 


The  techniques  and  results  on  the  laser  heating  strain- induced  (LHSI) 

, channel  waveguide  are  presented  in  Section  II.  The  automated  method  of 

^ fabricating  the  LHSI  waveguide  circuits  and  their  properties  are  given  in 

> Section  III.  The  results  of  introducing  the  residual  surface  strain  by 
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polishing  in  PLZT  and  LiNb03  and  the  techniques  for  measuring  the  effect  of 
strain  are  presented  in  Section  IV.  Section  V summarizes  the  progress  and 
results  made  in  the  course  of  this  study. 
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SECTION  II 

LASER  HEATING  STRAIN- INDUCED  CHANNEL 
WAVEGUIDES  IN  GLASS 


Introduction  of  residual  strain  in  a uniform  material  medium  can  cause  suffi- 
cient refractive  index  change  to  effect  optical  waveguiding.  To  explore  the 
utility  of  this  approach  for  fabricating  multimode- channel  waveguide  devices, 
we  have  developed  a refined  laser  heating  strain-induced  (LHSI)  technique 

using  filter  glass  wafers.  Our  approach  is  similar  to  that  reported  by  Pav- 
2 

lopoulos,  et  al  . First,  the  waveguide  dimension  in  our  case  is  of  the  order 
of  100/jm;  therefore  it  is  relatively  easy  to  couple  the  optical  beam  to  the 
waveguide.  More  importantly,  the  channel  waveguide  produced  by  the  present 
method  does  not  involve  burning  or  grooving  of  the  glass  wafer,  and  the  wave 
is  guided  in  a semicylindrical  channel  on  the  surface,  caused  by  the  localized 
refractive  index  gradient  due  to  the  laser  heating.  In  this  section,  the  tech- 
niques of  producing  this  type  of  channel  waveguide,  its  refractive  index  pro- 
file, and  waveguide  properties  are  presented. 


TECHNIQUES  OF  LASER  HEATING 


Properties  of  Glass  for  LHSI  Waveguides 


The  technique  of  laser  heating  and  machining  has  been  under  development 
since  the  invention  of  laser.  The  ability  of  achieving  high  energy  density, 
high  resolution  and  the  ease  of  beam  control  make  laser  heating  ideally 
suited  for  fabricating  miniature  or  microscopic  structures  such  as  encoun- 
tered in  the  microelectronics  industry.  To  extend  laser  heating  for  the  fab- 
rication of  multimode  integrated  optic  devices,  it  is  necessary  to  identify  a 
material  which  will  absorb  the  laser  energy  effectively  for  the  heating  laser 
beam  and  yet  transmit  the  propagating  laser  beam  with  low  loss.  We  have 


found  that  the  Corning  red  filter  glass  CS-2-60  satisfies  this  requirement. 
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The  properties  of  this  glass,  as  provided  by  Corning  Glass  are: 


Thermal 


Corning 

Expansion 

01<iss 

Glass 

Coefficient 

Number 

Type 

Dopant 

(x  10”  /cm°C) 

2408 

Zinc 

Soda 

Cd-Suifo- 

Selenade 

43 

Annealing 

St  iitenitij.! 

\S  orking 

Strain 

Temp.  f°C) 

rcniii. 

Temp.  (®C  ) 

Point  ( C) 

537 

7<i'l 

108  3 

501 

The  optical  properties  of  this  glass  were  measured  at  the  Argon  laser  and 
HeNe  laser  wavelengths.  For  the  Argon  laser  of  all  line  output  of  a Coherent 
Radiation  model  52G  laser  the  absorption  coefficient  was  measured  to  be 
67  cm“^.  The  refractive  index  and  absorption  coefficient  of  this  glass  at  the 

» i 

6328  A HeNe  wavelength  is  1.  507  and  0.  3 cm  respectively.  We  have  also 
measured  the  softening  temperature  by  differential  thermal  analysis  and  found 
that  it  is  betv/een  760-790°C.  The  wafer  is  polished  optically  flat,  with  a 
thickness  of  t5q)ically  2- 3mm,  and  overall  dimension  of  2.  5 x 5cm. 


Laser  Heating  and  Temperature  Profile 

To  produce  the  channel  waveguide,  the  sample  is  placed  near  the  focal  point 
of  an  Argon  laser  beam  using  a f/2  or  f/3  optical  lens,  and  is  mounted  on  a 
x-y  micropositioning  stage  driven  by  stepping  motors.  By  controlling  the 
drive  speed  and  laser  power  level,  channel  waveguides  of  various  dimensions 
were  obtained. 

The  temperature  rise  due  to  the  laser  heating  has  been  calculated  following 

4 

Pittaway's  treatment  . 

The  temperature  rise  T(”K)  in  a material  of  thermal  conductivity  K (watt/ 

2 

cm  ®C/cm)  due  to  a scanning  laser  beam  of  scanning  velocity  v (cm/sec) 
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focused  on  the  surface  of  a semi-infinite  target  with  a power  of  Eq  watts, 
and  a Gaussian  profile  of  radius  d cm  is  given  by: 


TT  / ^ 

j exp[{(m  - + B^]  sin^p  - ^^2^  (m^-3^)]dp  (1) 


where  Pd  and  md  are  the  distances  measured  from  the  center  of  the  beam,  per- 
pendicular and  along  the  scanning  direction  of  the  laser  beam  respectively, 
y = 4K /vd  and  K is  the  thermal  diffusivity  (cm^/sec). 


For  the  case  where  v = 0,  the  peak  temperature  is; 


T = E^/2Kd 


A 2 ~ 2 

using  typical  values  for  glass  (K  = 58  x 10*'*cm  /sec  and  K = 1.  18  x 10 
2 

watt/cm  “C/cm),  for  a laser  beam  of  available  power  at  the  medium  surface 
of  0.  5 watt  (value  used  for  our  experiments)  and  radius  of  d = 100pm,  we 
have 


T = 0.  5/1.  18  X 10"  • 2 \/n  = 1200“C 

P 


This  is  the  expected  peak  temperature  at  the  center  of  the  laser  scanned 
track.  Since  the  thermal  diffusion  time  for  a 100pm  radius  heated  spot  is 
approximately  given  by 


4x58x  10' 


= 4 msec. 


any  scanning  velocity  slower  than 


V ~ ^ ~ 2cm/sec 


can  be  considered  essentially  stationary. 

Calculation  of  temperature  rise  as  a function  of  m and  3 under  slow  scan 
condition  (y  ^ = 0)  was  performed  by  numerical  integration  and  the  result  is 
shown  as  a 3-d  plot  in  Figure  1.  The  normalized  temperature  at  3=0  is  also 
plotted  as  a function  of  m in  Figure  2.  It  is  seen  that  the  temperature- rise 
falls  to  about  50  percent  of  the  peak  value  at  m=l  or  at  a Gaussian  radius 
from  the  center  of  the  beam.  For  the  type  of  glass  used,  we  will  have  a track 
of  approximately  100/jm  radius  heated  above  the  softening  temperature. 


Surface  Deformation  and  Laser  Beam 

The  technique  of  laser  scanning  will  be  described  in  detail  in  a later  section. 
Basically,  the  optically  polished  glass  wafer  is  transported  across  the  focal 
point  of  an  Ar  laser  beam  and  a heated  track  is  produced.  Since  the  glass 
temperature  around  the  center  of  the  track  has  been  heated  above  the  soften- 
ing temperature  as  was  calculated  above,  and  due  to  thermal  expansion, 
this  heated  track  is  physically  raised.  Upon  subsequent  cooling,  a 
permanent  raised  track  results.  Figure  3(a)  shows  a typical  profile  of  such  a 
laser-heating-produced  ridge  as  was  measured  by  a Tallysurf  profilrecorder. 
A scanning  electron  micrograph  of  the  track  is  shown  in  Figure  3(b).  Since 
this  track  (or  ridge)  dimension  is  related  to  its  waveguiding  property,  we 
have  measured  the  ridge  height  h and  width  at  point,  w,  as  a function  of  the 
scanning  speed  and  laser  output  power  P,  (The  absorpted  power  is  approxi- 
mately 50  percent  of  the  output  power.  ) The  result  is  shown  in  Figures  4 and 
5.  The  scanning  speed  varies  from  4pm/sec  to  87pm/sec  and  the  Ar  laser 
all  line  output  power  varies  from  1.  04  to  0.  83  watt.  In  this  scanning  speed 
region,  the  temperature  rise  is  expected  to  be  essentially  independent  of  the 
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Figure  1.  Normalized  instantaneous  temperature  on  the 

glass  sample  by  a slow  scanning  focus  Ar  laser 
beam.  The  grid  spacings  are  half  of  the 
Gaussian  beam  radius.  Scanning  direction  is 
pointing  towards  lower  right. 
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Figure  2.  Normalized  temperature  as  a function  of  the 
normalized  distance  along  the  center  of  the 
scanning  track. 
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Profile  of  a laser  heated  ridge  on  Corning  filter  glass  CS-2-60 
(a)  Tallysurf  profile,  vertical  scale  0.  l/jm/ small  division, 
horizontal  scale  1 OO/um/major  division;  (b)  scanning  electron 
microscope  photograph  of  the  ridge,  magnification  500X. 


LASER  HEATED 


speed.  This  is  clearly  shown  in  Figure  5.  However,  the  detailed  shape  of 
the  heated  track  is  not  entirely  independent  of  the  scanning  speed,  as  is  seen 
in  Figure  4,  Since  temperature  rise,  therefore  the  product  h*  w=A,  is  propor- 
tional to  the  laser  power,  we  expect  A to  be  linearly  related  to  P.  This  is 
observed  in  Figure  5.  The  least-squares  fit  of  the  data  gives 

A = 320(P-0.  55>Mm^  (6) 

The  effect  of  laser  heating  to  the  localized  refractive  index  change  and  the 
waveguiding  properties  are  discussed  in  the  following  subsection. 


REFRACTIVE  INDEX  PROFILE 

Based  on  the  surface  deformation,  we  can  estimate  the  upper  limit  of  the 
elastic  strains  induced  in  the  neighborhood  outside  the  laser-heated  region, 
assuming  there  is  no  stress  relief  during  cooling  and  the  change  of  volume 
occurs  elastically.  For  h = l|um,  we  obtain  a strain  of  0.  005.  The  strain- 
induced  refractive  index  change  An^  at  this  point  can  in  turn  be  estimated  to 
be 

Ans  = - n^sp  ~ -0.  002  (7) 

where  n = 1.  507  is  the  measured  refractive  index  of  the  glass,  p is  the 
photoelastic  coefficient  and  is  measured  to  be  0.  2 (see  Section  IV).  Out- 
side this  region,  the  strain  is  tensile  and  falls  off  as  l/r,  as  does  An^. 

The  refractive  index  profile  can  be  understood  by  invoking  the  concept  of 

5 

the  fictive  temperature,  t,  and  the  relation  of  Tto  density  . In  the  central 
portion  of  the  laser-heated  region,  the  temperature  rise  exceeds  t^,  the 
existing  fictive  temperature  of  the  glass.  The  rate  of  cooling  and  consequently 
the  fictive  temperature  is  expected  to  be  dependent  on  the  radial  coordinate. 
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Since  the  center  of  the  heated  region  cools  at  a lower  rate  than  the  outer  i 

region,  it  will  acquire  the  lowest  rvalue,  highest  density  and  therefore  | 

the  highest  refractive  index.  Based  on  this  deduction,  the  cross-sectional  j 

profile  of  the  index  of  refraction  of  the  laser-heating- induced  waveguide  | 

should  take  a W shape.  The  central  hump  corresponds  to  the  region  of  | 

higher  density  resulting  from  a lower  r,  and  the  wings  are  due  to  the  strain-  | 

induced  effect. 

Experimental  measurements  of  the  refractive  index  profile  were  obtained 
using  a Jamin-Lebedeff  interference  microscope.  The  laser- written  track 
as  viewed  from  the  top  is  shown  in  Figure  6(a).  The  glass  slab  containing 
a cross  section  of  the  laser-heating- induced  waveguide  was  cut  and  polished 
optically  flat  and  parallel  on  both  sides  to  expose  the  cross  section.  A 
photograph  of  a typical  interference  pattern  of  the  waveguide  cross  section 
is  shown  in  Figure  6(b).  The  semicylindrical  symmetry  of  the  profile  is  evi-  j5 

dent  here.  Careful  examination  of  this  interference  pattern  as  the  analyzer  U 

II 

angle  on  the  microscope  is  varied  reveals  that  there  is  a refractive  index 
minimum  at  a distance  of  80pm  from  center.  The  final  analyzed  profile  for 
An  is  shown  in  Figure  7.  The  magnitude  of  the  refractive  index  change  in 

_ 3 

the  strained  region  is  about  3x10  , in  general  agreement  with  Eq.  (7). 

_ 3 

The  central  region  has  a localized  refractive  index  increase  of  An~  10 
Based  on  these  results,  we  can  estimate  the  numerical  aperture  N.  A.  of  the 
waveguide.  Using  the  analysis  for  step-index  optical  fiber  , N.  A.  = /2nAn, 
and  for  n = 1.  5,  An~  10  N.  A.  0.  05. 

The  high  residual  stress  could  be  substantially  reduced  by  a stress- relieve 
annealing  at  550°C  for  1 hour.  The  resultant  refractive  index  profile  in  the 
central  core  region  is  essentially  unaltered.  We  have  measured  the  incre- 
mental refractive  index  near  the  track  center  for  various  values  of  power 
P,  scanning  speed  v and  defocusing  condition.  The  results  are  shown  in 
Figures  8,  9,  and  10.  These  profiles  were  measured  along  a 45-degree  line 
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Interference  microscope  photograph  of  the  IJISI  track:  (a) 
top  view  where  the  fringes  are  mainly  due  to  the  raised  ridge 
(b)  cross-sectional  view  which  exposes  the  refractive  index 


DEFOCUSING  DISTANCE  (uM) 


Figure  10 


. Dependence  of  the  peak  value  of  the  refractive 
index  as  a function  of  the  defocussing  distance 
(distance  between  sample  surface  to  the  laser 
focal  point).  Data  obtained  for  100%  laser 
power  at  47pm /sec  scanning  speed. 
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from  the  center  of  each  track.  It  is  seen  that  by  proper  choice  of  P,  v and 
defocus ing  condition,  the  peak  value  of  An  can  be  controlled.  The  profile 
for  all  the  scanning  speed,  power  and  defocusing  condition  are  approximately 
parabolic.  The  waveguide  mode  analysis  to  follow  will  therefore  use  a para- 
bolic profile  approximation. 


RESIDUAL  STRESSES  IN  UISI  WAVEGUIDES 

In  this  subsection  the  residual  stresses  introduced  into  the  glass  by  laser 
heating  are  discussed.  As  was  discussed  earlier,  the  stresses  are  substan- 
tial and,  under  some  conditions,  can  cause  an  intriguing  form  of  fracture. 

It  should  be  noted  that  the  residual  stresses  cause  a reduction  of  the  refrac- 
tive index  and  acted  to  trap  the  guided  wave  which  propagates  in  the  softened 
region.  The  residual  stresses  can  be  substantially  removed  by  a stress 
relief  anneal,  the  laser  heated  region  still  acts  as  a waveguide  and  the 
material  is  no  longer  susceptible  to  fracture.  All  LUSI  waveguides  and  devices 
used  for  waveguiding  measurements  reported  here  were  annealed  and  are 
stable  and  fracture  free. 

The  residual  stresses  can  be  estimated  from  the  refractive  index  profile 
shown  in  Figure  7.  The  change  in  refractive  index  is  given  by  Eq.  (7). 

From  this  we  can  derive 

a = - (8) 

n'^ps 


where  E is  Young's  modulus.  Residual  radial  stresses  calculated  from  the 
index  measurements  made  on  the  8 0pm  radius  laser  heated  track  shown  in 
Figure  7 are  shown  in  Figure  11.  For  these  calculations  the  values  of  the 
constants  were  assumed  to  be  E = 9.  9 x 10  psi,  n = 1.  5 and  p ~ 0.  5.  The 
stresses  are  plotted  as  a function  of  the  distance  from  the  outside  of  the 
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heated  zone  normalized  by  a,  the  radius  of  the  heated  zone.  Note  that  the 
radial  stress  field  is  tensile  and  decreases  approximately  linearly  with  dis- 
tance from  the  heated  region. 

7 

These  stresses  can  be  estimated  with  the  aid  of  the  following  simple  model  . 
Assume  that  a cylindrical  region  of  radius  a,  in  the  glass  is  suddenly  heated 
to  above  the  softening  point,  (the  fictive  temperature).  Above  the 
cylindrical  region  is  stress-free.  As  the  region  cools  through  t ^ it  becomes 
an  elastic  solid  and  starts  exerting  a radial  tensile  stress  on  the  surrounding 
material  as  it  shrinks.  The  radial  tensile  stress  field  resulting  from  the 
cylinder  cooling  to  room  temperature  is  given  by 


a EAT  ,a 
r = ■iri-vT 


(9) 


where  a is  the  coefficient  of  thermal  expansion,  v is  Poisson's  ratio  and  r is 

0 

the  radial  distance  from  the  axis  of  the  cylinder;  E = 9.  9 x 10  psi,  v - . 2, 

“6  “1 

and  O'  is  assumed  to  be  6 x lO"  °C  . The  radial  stresses  calculated  from 
the  simple  shrink  fit  model  are  shown  in  the  solid  line  on  Figure  11.  The 
model  does  not  predict  the  observed  linear  dependence  of  with  r but  the 
values  of  stress  are  of  the  right  order.  It  should  be  pointed  out,  however, 
that  the  data  on  both  the  Young's  modulus,  K,  and  tlie  thermal  expan- 
sion coefficient,  at  are  estimated,  since  values  measured  at  room  tempera- 
ture are  not  expected  to  be  valid.  For  this  reason  the  magnitude  of  stress 
plotted  on  Figure  11  should  only  be  considered  as  approximate. 

In  the  above  discussion  we  have  only  considered  the  stresses  in  the  material 
outside  the  heated  region.  The  stresses  inside  the  heated  region  cannot  be 
determined  from  interference  microscope  measurements  such  as  those 
shown  in  Figure  7 because  heating  above  the  softening  point  has  also  altered 
the  density.  In  this  case  the  index  of  refraction  is  a function  of  both  the 
density  and  the  residual  stress.  The  two  cannot  be  separated. 
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When  an  unannealed  waveguide  comes  into  contact  with  a foreign  object  (e.  g.  , 

a coupling  prism),  it  sometimes  fractures  in  a most  unusual  manner.  A 

crack  forms  at  the  point  of  contact,  usually  at  the  edge  of  the  glass  where  ’ 

the  waveguide  terminates,  and  runs  along  one  side  of  the  guide  parallel  to  it. 

After  propagation  has  continued  for  some  distance  the  crack  turns  and  runs  ’ 

i; 

through  the  guide  at  right  angles  to  it  until  it  reaches  the  opposite  side.  At  i 

this  point  the  crack  momentarily  stops  and  then  starts  moving  again  parallel  | 

to  the  guide  as  before  but  now  on  the  opposite  side.  After  running  approxi-  | 

mately  the  same  distance  as  before  the  crack  turns  and  runs  across  the  guide  | 

and  stops.  At  this  point  the  process  repeats  itself.  A top  view  of  this  frac-  | 

ture  process  is  shown  in  Figure  12.  The  crack  itself  is  cylindrical,  as  ex-  { 

pected  from  the  radial  tensile  stress  field  shown  in  the  interference  measure- 
ments.  A scanning  electron  micrograph  of  the  end  view  of  the  crack  is  shown  j 

in  Figure  13.  The  ridge  running  down  the  surface  of  the  glass  is  the  top  of  i 

the  region  heated  above  the  softening  point  by  the  moving  argon  laser  beam. 

The  crack  runs  at  quite  a slow  rate  and  it  takes  many  minutes  to  traverse  the 
whole  length  of  the  guide.  Thus,  the  fracture  is  subcritical  and  occurs 
most  likely  by  a stress  corrosion  mechanism. 

In  what  follows,  a phenomenological  model  describing  the  fracture  is  pre- 
sented. Quantitative  discussion  must  postponed  till  further  analysis  is 
made.  As  we  have  seen,  the  waveguide  is  a semicylindrical  region  at  the 
glass  surface  that  was  heated  above  the  glass  transition  temperature  by  a 
moving  argon  laser  beam.  The  softened  region  is  rapidly  quenched  following 
passage  of  the  laser.  The  quenched  region  has  a lower  density  than  the  bulk 
glass  as  is  evident  from  the  ridge  seen  in  Figure  13.  Heat  flow  during  cool-  ; 

ing  is  predominantly  into  the  glass  and  the  resulting  residual  stress  field  has  ' 

radial  symmetry  about  the  axis  of  the  guide. 


The  principal  stresses  in  this  case  are  a^,  a radial  tensile  stress, 
a compressive  tangential  ijAess,  and  a tensile  stress  lying  along  the  axis 
of  the  guide.  The  principal  stresses  are  shown  on  Figure  14  which  is  a 
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Figure  14.  Schematic  diagram  of  residual  stresses  showing  model 

of  fracture  experienced  by  the  unannealed  LHSI  waveguide. 


schematic  of  the  proposed  fracture  model.  The  magnitudes  of  the  principal 
tensile  stresses,  and  depend  on  the  quenching  rate.  Glass  is  a visco- 
elastic material  and  can  undergo  stress  relaxation  during  cooling,  de- 
pends on  heat  flow  along  the  guide  in  the  z direction  which  in  turn  depends 
on  the  rate  of  traverse  of  the  laser  beam.  We  presume  that  a >a  as  indi- 
cated  on  the  figure.  The  initial  crack  will  be  perpendicular  to  and  will  run 
parallel  to  the  guide  in  the  z-direction.  is  a residual  stress  and  is 
relieved  as  the  crack  propagates.  Since  the  crack  will  follow  the  direction  of 
maximum  tensile  stress  it  will  turn  across  the  guide  as  the  ratio 
decreases.  When  the  crack  runs  across  the  guide  is  relieved  and  the  crack 
halts.  We  have  observed  that  the  crack  bifurcates  (i.  e.  , branches)  as  it 
stops.  One  of  the  branches  is  in  a favorable  orientation  with  respect  to 
in  the  uncracked  region  and  the  process  repeats  itself. 

A number  of  questions  concerning  the  model  remain  to  be  answered.  First, 
as  mentioned  above,  due  to  uncertainties  in  the  Young's  modulus  and 
in  the  thermal  expansion  coefficient,  the  stresses  in  the  region  outside  the 
heated  zone  where  the  crack  is  running  are  not  known  with  much  certainty. 
Another  estimate  of  the  radial  stress  can  be  obtained  from  the  separation 
at  the  crack  (i.  e.  , at  point  A in  Figure  13).  The  separation  was  measured 
by  scanning  electron  microscopy  to  be  about  Ip  for  the  crack  running  235pm 
outside  of  the  center  of  the  guide.  This  corresponds  to  a surface  radial  strain 

_ O 

of  2.  1 X 10  . The  residual  stress  corresponding  to  this  strain  is  21000  psi 

which  is  on  the  order  of  the  stresses  shown  in  Figure  11. 

While  the  magnitude  of  the  residual  stress  is  reasonable  the  position  of  the 
crack  is  not.  The  maximum  radial  tensile  stress  is  at  the  interface  between 
the  heated  zone  and  the  surrounding  material.  The  observed  fractures  occur 
at  a considerably  distance  from  the  heated  zone  as  obvious  from  Figure  13. 
Typically  the  distance  of  the  crack  to  the  waveguide  is  an  order  of  magnitude 
larger  than  the  radial  dimensions  of  the  guide.  At  this  distance  from  the 
guide  the  stress  birefrigence  shown  in  Figure  7 is  small  indicating  a low 
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stress  level.  At  this  time  we  do  not  understand  why  the  guides  crack  where 
they  do. 


Alternate  models  of  the  residual  stress  in  the  glass  following  passage  of  the 
laser  beam  have  been  suggested  which  more  realistically  treat  the  glass  as 
a viscoelastic  material®.  More  analysis  and  experimentation  are  needed  to 
achieve  a better  understanding  of  this  crack  mechanism.  It  should  be  empha- 
sized here  that  in  all  the  optical  waveguiding  experiments,  the  LHSI  wave- 
guides are  all  preannealed  and  stresses  are  substantially  removed.  No 
crack  was  observed  in  the  annealed  samples. 

LHSI  WAVEGUIDE  PROPERTIES 
Waveguide  Mode  Theory  for  LHSI  Guides 

The  LHSI  guides  described  above  may  be  conveniently  modeled  as  having  an 
index  profile  with  semicircular  symmetry.  That  is,  there  is  a line  of 
symmetry  on  the  surface  of  the  guide  material  with  the  index  of  the  guide  a 
function  only  of  distance  from  this  line  inside  the  material.  In  this  section 
we  present  a simple  analysis  of  the  modes  of  multimode  graded-index  fibers 
(circular  symmetry)  and  semicircular  graded-index  guides.  Normalized 
parameters  are  defined  so  that  a single  equation  or  a single  universal  chart 
describes  all  such  guides  having  a particular  index  profile  shape,  such  as  para- 
bolic, Gaussian,  etc.  The  initial  steps  in  the  problem  description  come  from 
Reference  6.  We  first  consider  circularly  symmetric  graded-index  fibers 
and  will  later  explain  the  modifications  necessary  for  a description  of  semi- 
circular guides. 

We  assume  that  the  fiber  waveguide  is  circularly  symmetric  with  a small 
index  change  in  the  core  region,  and  that  index  variations  on  the  scale  of  a 
wavelength  are  small.  The  index  profile  is  a continuous  function  of  radial 


position,  n(r).  At  any  radial  position  r the  wave  vector  has  the  magnitude 
k(r)  = 2nn(r)/A.  = k^n(r). 


The  modes  of  this  waveguide  are  described  by  two  integer  mode  numbers,  v 
and  u.  The  azimuthal  mode  number,  v describes  the  number  of  cycles  of  the 
electric  field  in  the  azimuthal  direction.  At  position  r the  azimuthal  compon- 
ent of  the  wave  vector  is  therefore  v/r  ( = 0,  1,  . . . ) in  order  that  the  phase 
go  through  V cycles  around  the  circumference  2Tir;  e.  g.  , the  variation  of 
amplitude  with  azimuthal  angle  0 is  either  sin  (v0/r)  or  cos  (vS/r).  If  a 
mode  exists  which  propagates  down  the  fiber  with  wave  vector  0,  then  by  a 
geometrical  construction  the  radial  component  of  the  wave  vector  at  radial 
position  r is 


u(r) 


r k^(r) 


(10) 


Only  between  two  radial  positions  Rj  and  R2  is  u(r)  real  and  the  radial  phase 
variation  periodic.  u(R^)  = u(R2)  = 0 and  R^  < R2.  Outside  this  region  the 
field  is  aperiodic.  The  wave  can  be  through  to  undergo  reflections  from  the 
graded  index  profile  at  R^  and  R^.  Such  a reflection  introduces  a phase 
shift^  of  tt/4.  The  total  radial  phase  change  between  R^  and  R2.  including 
these  phase  shifts,  must  add  up  to  an  integer  multiple  of  tt  for  a mode  to 
exist,  leading  to  the  mode  dispersion  equation 


9 1/2  . 

^ - B^]  dr  = (p  (p  = 0,  1,  2,  . . . ) 

r2  2 


(11) 


A similar  mode  dispersion  equation  applies 


for  diffused  planar  guides 
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Let  the  index  profile  be  written  as 


n(r)^  = nj^^  + (n^^  - n^j^)f(r/R^)  (12) 

where  n is  the  index  at  the  center  of  the  fiber,  n.  is  the  constant  index  away 

3.  ^ 

from  the  graded  core,  and  f(r/R^)  is  the  normalized  shape  of  the  index  pro- 
file. f(0)  = 1 and  f(r/R  ) declines  monotonically  to  zero  with  increasing  r. 

° 2 
For  example,  a parabolic  index  fiber  would  have  f(u)  = 1 - u for  0 < u < 1. 

R is  the  scale  factor  for  the  index  shape, 
o 

Two  normalized  parameters  can  be  defined:  The  normalized  guide  radius  V, 
and  the  normalized  mode  effective  index  b 


V = k R [n  ^ - n. 

o o a b 


n i-r  "n, 

eff  b 

2 2 
n -n, 
a b 


where  the  effective  index  of  a waveguide  mode  is  n^^^  = 6/k^.  Similar  nor- 
malized parameters  can  be  used  to  define  thin  film  and  diffused  planar  wave- 
9,10 

guides 

The  mode  dispersion  Eq.  (11)  cm  now  be  written  in  normalized  form  as 


,2  1/2 


f rv^f(u)  - V^b  - — ] du  = (p  + 1/2)tt 
J 1 1 2 


where  we  have  made  the  substitution  u = r/R^,  and  the  limits  and  are 
where  the  integrand  equals  zero,  with  < U2.  For  a given  profile  shape  f(u), 
this  normalized  mode  dispersion  Eq.  (15)  describes  the  modes  of  all  fibers 
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having  such  a shape  in  terms  of  the  normalized  parameters  b and  V.  Equiv- 
alently, for  each  shape  f(u)  a universal  chart  of  b versus  V can  be  plotted 
from  Eq.  (15)  which  describes  the  modes  of  all  waveguides  having  this  index 
shape. 

The  total  number  of  modes  in  a graded- index  fiber  is  shown  in  Reference  6 to 
be 


\ 


i 


M 


rn(r)^ 


nj^^Jrdr 


(16) 


where  a is  where  n(r)  reaches  the  cladding  index  n^,  and  M includes  both  a 
two-fold  degeneracy  in  polarization  and  a two-fold  degeneracy  in  azimuthal 
orientation.  Changing  variables  to  u = i’/Rq  and  substituting  Eqs.  (12)  and 
(13)  gives 


M 


V‘ 


f 

> o 


f(u)udu 


(17) 


Many  of  the  index  profiles  for  ITISI  guides  can  be  approximated  as  parabolic. 
We  now  show  how  a parabolic  index  profile  leads  to  an  extremely  simple  form 
for  the  mode  dispersion  equation. 

For  a parabolic  index  fiber,  the  shape  function  f(u)  is 
I"  1 - u^  for  0 < u < 1 

f(u)  (18) 

; 0 for  u > 1 

Substituting  this  into  the  mode  dispersion  Eq.  (15),  that  equation  can  then  be 
integrated.  The  mode  dispersion  equation  for  a parabolic  index  fiber  becomes 
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V(1  -b)  = 2(2/u  + v+l)  = 2m  (m  = 1,  2,  . 


(19) 


. . ) 

where  all  modes  with  indices  ^ and  v giving  the  same  value  of  m = 2/u  + v + 1 
are  degenerate.  This  equation  can  be  shown  to  be  equivalent  to  the  unnormal- 
ized mode  dispersion  equation  derived  in  Reference  11  for  parabolic  index 
fibers.  A universal  chart  of  b versus  V describing  the  modes  of  all  para- 
bolic index  fibers  can  easily  be  plotted  and  is  shown  as  Figure  15. 

For  a parabolic  index  fiber  the  total  number  of  modes  M can  be  written  in 
normalized  form  using  (17)  as: 


M includes  both  the  two-fold  degeneracy  in  polarization  (x  or  y)  and  the  two- 
fold degeneracy  in  orientation  (azimuthal  variation  of  sin  (v6/r)  or  cos 
(v9/r)). 

The  above  results  may  be  used  to  describe  parabolic  index  profile  fibers. 
These  fibers  are  of  great  importance  because  of  the  small  variations  in  group 
velocity  for  the  various  modes,  leading  to  a narrow  impulse  response  and 
the  capability  for  transmitting  information  at  high  data  rates. 

To  apply  the  results  given  by  Eq.  (15)  and  (17)  for  the  general  graded-index 
fiber  to  the  semicircular  model  of  LHSI  guides,  we  note  that  a semicircular 
guide  is  one  which  looks  like  half  of  a circularly  symmetric  fiber  which  has 
been  sliced  along  its  axis.  Denoting  the  plane  of  this  cut  as  y = 0 and  assum- 
ing the  guide  is  bounded  above  by  air,  the  index  profile  in  rectangular 
coordinates  is: 

o ^ f((x^  + y^)/D^)  for  y^  0 

n(x,y)^  ^ b (21) 

I 1 for  y < 0 
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In  cylindrical  coordinates,  the  index  profile  is  given  by  (12)  and  (21)  when 
the  azimuthal  angle  is  in  the  range  0 < 6 < n.  Outside  the  range  of  6,  n(r)  = 1. 

At  the  y = 0 surface  of  this  semi-circular  guide,  the  electric  field  amplitude 
of  any  mode  is  small  compared  to  its  maximum  inside  the  guide,  but  a weak 
evanescent  field  still  exists  for  y < 0.  The  larger  n^  is,  the  weaker  this  field. 
If  we  make  the  assumption  that  the  y = 0 surface  is  a perfect  reflector  (or 
that  n^  is  large,  or  that  n = 0 for  y < 0),  then  the  field  amplitude  of  each  mode 
is  zero  across  the  y = 0 plane,  and  we  can  make  use  of  the  previous 
description  of  the  modes  of  a circular  fiber.  Modes  of  the  fiber  having  a null 
along  y = 0 thus  are  modes  of  the  semi-circular  guide.  Modes  with  v = 0 
are  ruled  out,  but  all  modes  with  other  values  of  v have  null  planes.  The 
requirement  that  this  null  plane  be  along  y = 0 merely  serves  to  remove  the 
degeneracy  in  azimuthal  orientation  previously  described.  Only  modes  with 
azimuthly  variation  sin  (v  fi/r)  with  v > 1 are  permitted,  and  not  modes  with 
cos  (V  0/r)  variation  or  v = 0.  Equation  (15)  may  still  be  used  to  find  the 
allowed  combinations  of  b and  V for  a particular  f(u)  and  mode  indices 
v>  1,  /u  > 0. 

The  total  number  of  modes  in  a semi-circular  guide  is  halved  from  that 
given  by  Eq.  (12)  due  to  the  elimination  of  the  orientational  degeneracy,  and 
is  further  reduced  by  not  permitting  modes  with  v = 0. 

For  a parabolic  index  fiber,  Eq.  (20)  gives  the  total  number  of  modes,  in- 
cluding the  four-fold  degeneracy.  A number  of  possible  pairs  of  p and 
for  a particular  V is 

N = 1 M - -Sr  (22) 

4 In 

For  a semi-circular  guide,  we  subtract  the  number  of  pairs  with  v = 0.  For 
each  value  of  m there  is  one  such  pair,  so  the  number  of  pairs  eliminated 
equals  the  maximum  value  at  m which  is  V.  Multiplying  by  the  two-fold 
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degeneracy  in  polarization,  we  obtain  the  total  number  of  modes  for  a para- 
bolic semi-circular  guide  as 

M = ^ - 2V  (23) 

For  a typical  LHSI  guide,  = SOpm,  n^^  = 1.  51.  and  (n^  - n^^)  = 1.  0 x lO'^, 

so  V ~ 44  for  operation  at  X ^ = 6328A.  This  leads  to  a total  number  of  modes 
M T 151. 


Losses  Due  to  Index  Reversal  in  LHSI  Guides 


The  LHSI  optical  waveguides  have  a higher- index  central  region  with  the 
index  decreasing  with  distance  from  the  axis  of  symmetry.  However,  out- 
side a certain  radius  the  index  profile  reverses  and  begins  to  increase, 
reaching  a value  in  the  bulk  material  which  is  often  higher  than  in  the  central 
guiding  region.  Energy  in  the  waveguide  would  be  expected  to  tunnel  through 
the  low  index  region  and  leak  into  the  higher-index  bulk  material.  We  pre- 
sent here  calculations  which  show  that  the  losses  due  to  this  leakage  should 
be  very  small  for  the  LHSI  guides. 


Wo  again  model  the  guides  as  having  semi-circular  symmetry  and  a parabolic 
profile,  so  that  we  can  adapt  results  obtained  for  parabolic-index  fibers. 
Marcuse^  ^ shows  that  a parabolic-index  fiber  defined  by; 

n(r)^  = n^^  (1  - 2Ar^/a^)  (24) 

has  modes  with  transverse  electric  field  amplitudes  in  Cartesian  coordinates 
of 


E (x,  v)  = 
mn  ’ • 


AH  ( 
m 


\/2x/W)H  ( /2y/W)exp  l-(x' 


y^)/w 


2 


(25) 
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where  m and  n are  the  mode  indices,  A is  a constant,  denotes  a Hermite 
polynomial,  and  W can  be  considered  the  mode  radius  given  by 

W - (2/A)^^^  (a/n^k)^^^  (26 

These  modes  have  the  same  form  as  the  modes  of  a laser  resonator,  except 

that  in  the  latter  case  W is  a function  of  axial  position  as  well.  The  electric 

1 3 

fields  of  the  modes  can  also  be  written  in  cylindrical  coordinates  as  : 

E = (2r^/W^)exp(-r^/W“)  r q«]  (27 

pq  o p 

where  p and  q are  the  radial  and  azimuthal  mode  indices,  respectively,  E 

, q 

is  a constant,  and  l_,  is  the  generalized  Laguerre  polynomial. 


Using  the  typical  LHSI  profile  of  Figure  7,  we  find  it  is  approximately  para- 
bolic, and  with  n^  = 1.  51  we  calculate  that  W "=  17pm.  Thus  the  fields  at  the 
index  reversal  radius  of  about  = 8Qum  are  expected  to  be  small  and  the 
leakage  into  the  higher-index  bulk  material  should  be  very  slow. 

VVe  can  make  an  approximate  calculation  of  this  loss  rate  by  modifying  the 

14 

theory  of  Marcatili  and  Miller  , used  to  calculate  the  radiation  loss  from 
curved  guides.  They  approximate  the  attenuation  coefficients  as: 


where  2Z^  is  the  colimated-beam  length  associated  with  the  transverse  field 
distribution,  i.  e. , it  is  the  length  in  which  the  beam  radius  would  double  due 
to  diffraction  in  material  of  index  n^.  For  the  case  considered  here 


2Z  = n {2Wri\ 
c o 


I 


and  for  the  profile  of  Figure  7,  2Z^  '=  2.  8mm.  The  ratio  is  the  frac- 

tion of  the  mode  energy  in  the  region  where  loss  occurs.  For  loss  due  to 


index  reversal,  we  can  consider  all  energy  is  lost  for 
then  given  by: 


e j are 


” 2TT 

^ rdr  ( E(r,  ij)^d 

" R o 

o 


(30) 


2ri 


= J I'dr  J E(r,  gi)' 


(31) 


where  E(r,  $)  is  given  by  Eq.  (27). 

By  approximating  the  Laguerre  polynomial  in  Eq.  (27)  for  large  values  of 
x/W  and/or  small  values  of  q,  we  find  that  approximately: 


m 


T ~ 


-u 


m-i 

o Y ''o 

Z. 


(32) 


i=0 


where  u = 2R  and  m = 2p+q. 

o o ^ ^ 


We  have  calculated  a from  Eq.  (32)  using  W - 17pm,  n^  = 1.  51,  X = 0.  6328pm, 
and  Rq=  80pm  as  in  Figure  7.  The  results  are  shown  in  Table  1 for  values  of 
m from  0 to  28.  While  these  calculations  are  very  approximate,  they  do  show 
that  the  loss  rate  should  be  very  low  for  most  waveguide  modes.  This  is  ex- 
pected for  the  small  value  of  W compared  to  R . 


t. 
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Table  1.  Approximate  Values  of  Loss  Coefficient  a for 
Various  Mode  Indices  [m  = 2p  + l due  to  index 
reversal  in  the  guide  of  Figure  7 using  Eq.  (32)]. 


m = 2p+l 


Approximate  ty 
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10-20, 

2 

6 

X 
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1 

o 

4 

1 
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6 

7 
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8 

3 

X 
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10 
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22 
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24 
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X 
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26 

2 
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28 

6 
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Properties  of  LHSI  Waveguide 

The  waveguiding  properties  of  the  LHSI  glass  waveguide  has  been  measured 

using  a glass  prism  coupler  with  a base  angle  9 = 45°  and  a refractive  inde> 

P 

The  parameters  measured  are  the  waveguide  acceptance  angle 


of  n =1.  8. 

P 

and  the  waveguide  losses. 
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The  effective  index  of  refraction  of  the  waveguide  n^^^  is  related  to  the  prism 
input  coupling  angle  0^  by  the  Snell's  law. 

"eff  = + 

where  n sin  0„  = sin  0, 

p 2 1 

From  this,  we  have 


dn  „„  c^s0, 


P’or  0p  = 45°,  n^  = 1.8  and  in  the  region  where  0^  = 22°,  we  have 


0.  516 


The  measurement  of  the  acceptance  angle  is  indirectly  obtained  by  plotting 
the  laser  beam  intensity  from  an  output  prism  as  the  incident  angle  0^^  is 
varied.  The  result  for  both  the  TE  and  TM  modes  is  shown  in  Figure  16. 

The  acceptance  angle  at  half  intensity  point  for  both  cases  is  8 minutes  or 
0.  0023  radian.  Therefore 

An  »»  = 0.  0023  - 0.  516  = 0.  0012  ( 

eii 

1 6 

Using  the  relation  of  numerical  aperture  N.  A.  for  optical  fiber  with  a 
step  index  An  of  core  index  n. 


2nAn 


RELATIVE  INTENSITY  ( Arb.  Units) 


Figure  16.  Relative  beam  intensity  at  a fixed  point  in  the  LHSI  wave- 
guide as  a function  of  the  incident  angle  0^  to  a 45 -degree 

glass  prism  (n  = 1.8)  coupler  for  both  TM  and  TK  modes. 
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we  can  estimate  the  N.  A.  of  the  LHSI  waveguide.  We  obtain 

N.  A.  = y 2 • 1.  5 • 1.  2 X 10"^  = 0.  06  (38) 

The  Angff  values  calculated  from  the  acceptance  angle  is  in  agreement  with 
the  An  value  observed  from  the  interference  microscope  measurement  of 


It  is  important  to  note  that  the  LHSI  waveguide  should  be  relatively  easy  to 
couple  to  typical  multimode  optical  fiber  which  possess  a N.  A.  value  of  0.  1- 
0.  2 and  a fiber  diameter  comparable  to  the  LHSI  waveguide  dimension. 

The  attenuation  of  the  LHSI  waveguide  was  also  measured  using  a sliding  out- 
put prism  coupler.  The  measured  relative  beam  intensity  as  a function  of 
waveguide  length  for  both  the  TE  and  TM  waves  is  shown  in  Figure  17.  The 
data  were  least-squares  fitted  with  an  exponential  curve  and  the  result  is 
also  shown  in  the  same  figure.  From  these  curves,  we  obtain  the  optical 
losses  as  follows: 


Optical  losses  for  TM  waves  = 1.  66db/cm 
Optical  losses  for  TE  waves  = 3.  Idb/cm 

These  results  as  well  as  the  LHSI  refractive  index  profile  has  been  reported 

12 

m a recent  publication 

Based  on  the  understanding  of  the  fictive  temperature  change  for  the  forma- 
tions of  the  waveguiding  region,  we  could  expect  an  increased  n and  there- 
for N.  A.  if  the  glass  wafer  were  preannealed  for  a longer  time  duration. 
This  approach  was  experimentally  investigated.  We  learn  that  the  annealing 
schedule  for  the  CS-2-60  glass  as  prepared  at  Corning  was  as  follows:^ 

• Heating  from  room  temperature  to  annealing  temperature  within 
15-20  minutes. 
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RELATIVE  BEAM  IN7 
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LENGTH  OF  WAVEGUIDE  (cm) 


Figure  17.  Relative  beam  intensity  as  a function  of  the  LIISl  waveguide 
path  length  for  both  the  TM  and  TTO  modes. 
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• stay  at  this  temperature  for  about  1/2  or  1 hour. 


• Cool  back  to  strain  point  501°C  at  a rate  of  10“ /minute. 

• Cool  to  room  temperature  rapidly  with  furnace  off. 

To  attempt  to  decrease  the  fictive  temperature,  we  have  chosen  a slower 
cooling  rate.  The  schedule  adapted  was  as  follows; 

• Heating  from  room  temperature  to  600°C  in  1/2  hour. 

• Stay  at  this  temperature  for  2 hours. 

• Cool  back  to  500°C  at  6°/hour  or  0.  l“/minute. 

• Cool  to  room  temperature  with  furnace  off. 

Using  two  glass  wafers  cut  from  the  same  plate  from  Corning,  we  prepared 
one  according  to  our  new  preannealing  schedule  and  leaving  the  other 
untreated.  LHSl  waveguides  were  then  made  on  both  wafers  under  identical 
conditions.  The  measured  waveguide  acceptance  angle  for  the  preannealed 
sample  is  shown  in  Figure  18.  It  is  seen  that  the  half  power  aperture  for  the 
TE  and  TM  modes  are  9.  5 and  10.  5 minutes  respectively,  as  compared  to 
the  as-received  sample  which  gives  an  half  power  aperture  of  8.  5 minutes 
for  both  the  TE  and  TM  wave.  This  indicates  that  preannealing  is  effective 
for  increasing  the  numerical  aperture  of  the  LHSI  waveguide. 

This  result  is  a further  indication  that  the  fictive  temperature  of  the  glass 
plays  an  important  role  in  the  refractive  index  of  the  glass  as  was  previous- 
ly expected.  However,  our  result  shows  the  range  of  modification  in  An  by 
slow  preannealing  is  not  very  large. 
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SECTION  III 

FABKICATION  OF  I.HSI  WAVEGUIDE  CIKCLTT  ELEMENTS 


FABRICATION  TECHNIQUE 
The  Laser  Beam  Exposure  System 

The  use  of  a laser  beam  for  LHSI  waveguide  fabrication  offers  a unique 
capability  in  the  realization  of  complex  optical  waveguide  devices.  To  im- 
plement this  capability,  a laser  beam  exposure  system  (LBES)  was  designed, 
constructed  and  tested.  In  this  system  the  focused  laser  beam  is  stationary 
and  the  sample  receiving  the  exposure  is  mounted  on  an  x-y  micropositioning 
stage,  controlled  by  a microcomputer.  An  electronic  shutter  is  also  incor- 
porated to  control  the  beam  on-off  duration.  This  system  has  been  used  for 
the  fabrication  of  channeled  waveguide  bends,  mode  selectors,  and  couplers. 

It  should  be  applicable  for  multimode  optical  waveguide  device  fabrication 
where  the  waveguiding  region  is  produced  by  the  photon  energy  in  the  laser 
beam. 

The  schematic  diagram  of  the  LBES  is  shown  in  Figure  19.  The  laser 
employed  for  this  system  is  a Coherent  Radiation  Model  52G  Ar  laser  with 
an  maximum  all-line  power  output  of  about  1 watt.  The  laser  beam  is 
filtered  by  a spatial  filter,  controlled  by  an  electronic  shutter  and  is  focused 
by  a 30mm  f.  1.  lens  operating  at  an  F number  of  2 to  3.  The  x-y  micro- 
positioning stage  consists  of  two  Aerotech  model  ATS-302A1  linear  translators 
with  stepping  motor  drive,  driven  by  Aerotech  model  200  drivers.  A 
Motorola  M6800  microprocessor  (uP)  was  used  to  control  the  drivers  of  the 
microposition  stages.  The  schematic  diagram  of  the  electronic  circuits  is 
shown  in  Figure  20,  where  the  mnemonics  are  given  in  Table  2, 
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Table  2.  Signal  Mnemonics 


Y -Stage 

X -Stage 

PAO  - Home  Selected  (STAT) 

PHO  - Home  Selected  (STAT) 

PA  1 - Acknowledge  Home  Selected 

PHI  - Acknowledge  Home  Selected 

(CMD) 

(CMD) 

PA2  - Limit  Switch  Active  (STAT) 

PH2  - Limit  Switch  Active  (STAT) 

PAS  - Drive  FWD/H.EV  (CMD) 

PBS  - Drive  FWD/KEV  (CMD) 

PA4  - Direction  of  Travel  (STAT) 

PH4  - Direction  of  Travel  (STAT) 

PH5  - Laser  Shutter  Control  (CMD) 

PA5  - Spare 

PA6  - Spare 

PB6  - Spare 

PA7  - Spare 

PB7  - Spare 

CAl  - Spare 

CHI  - Spare 

CA2  - Step  Strobe  (\’ -Stage) 

CH2  - Step  Strobe  (X -Stage) 

In  the  following,  the  hardware  components  and  their  functions,  as  well  as  the 
software  and  programs  utilized  are  discussed  in  detail. 

The  drives  consist  of  two  independent  electronic  packages  capable  of  opera- 
tion in  a lical  or  remote  mode.  When  operated  in  the  local  mode,  the 
direction,  slew,  on/off  and  a single  step  function  are  controlled  by  switches 
on  the  front  panel  of  the  drives.  A potentiometer  adjusts  the  internally 
generated  stepping  rate  from  a minimum  of  55  steps/sec  (27.  5um/sec)  to  a 
maximum  of  500  steps  (250um/sec),  An  L.EU  lights  when  the  range  limit 
switch  has  been  actuated  indicating  the  stage  limit  of  travel  has  been  reached. 
The  range  limit  switch  also  removes  stepping  motor  drive  and  the  stage 
stops.  When  the  stepping  rate  was  varied,  it  was  determined  that  one  stage 

• resonated  at  192  steps/sec,  110  steps/sec,  and  71  steps/sec.  Measurements 

of  the  resonant  frequencies  were  not  made  on  the  second  stage  { \ axis) 
although  similar  resonant  points  were  evident. 
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When  the  stages  are  switched  to  tiie  remote  mode  of  operation,  the  M6800 
microprocessor  (uP)  provides  a 7usec  stepping  trigger  pulse  to  the  stage 
driver.  The  stepping  pulse  triggers  a one-shot  multivibrator  wliich  stretches 
the  pulse  to  the  30usec  width  required  by  the  driver  circuitry.  Th  stepping 
rate  in  the  remote  mode  is  entirely  dependent  on  the  rate  at  which  the  uP  ^ 

supplies  the  trigger  pulses.  The  maximum  pulse  rate  in  the  remote  mode  is  - 

approximately  the  same  as  in  local  mode,  although  the  aP  is  capable  of  i 

supplying  step  pulses  at  a rate  much  higher  than  the  stepper  motor  is  capable  1 

of  stepping.  j 

j 

•1 

1 

The  aP  interface  also  provides  the  exchange  of  status  signals  which  will  j 

indicate  whether  the  "home  select"  (initializes  origin  point)  switch  has  been 

actuated,  the  travel  limit  switches  have  been  actuated,  and  the  direction  the 

stage  is  stepping.  Command  signals  sent  to  the  driver  stage  consist  of  the 

step  strobe,  the  "acknowledge  home  selected"  command  which  will  reset  the 

home  select  F/F,  the  drive  forward/ reverse  signal,  and  a command  which 

will  open  or  close  the  laser  shutter. 

The  microcomputer  consists  of  an  M6800  microprocessor  situated  on  a - 

Motorola  MKK6800D1  evaluation  PCH  along  with  a ROM  (MCM6830L7, 

MIKHUG)  and  associated  HAM  (MCM6810),  a peripheral  interface  adapter 
(PIA,  MCM6820)  tor  communication  with  a Teletype  (TTY)  or  other  RS232 
compatible  communication  device,  and  a I’lA  for  communication  with  the 
translation  stages. 

A separate  memory  card  was  constructed  which  expanded  the  basic  evaluation 
PCH  memory  capacity  from  640  words  to  1536  words.  This  card  is  isolated 
from  the  evaulation  F^CR  bus  by  a bus  expander  card  (using  National  8833 
party  line  transceivers)  which  will  allow  up  to  100  MOS  devices  to  be  placed 
on  the  buffered  bus. 
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The  PIA  used  for  translation  stage  interfacing  has  16  parallel  programmable 
I/O  lines  and  4 control  lines.  Two  control  lines  are  strictly  input  and  can  be 
programmed  by  the  uP  to  generate  an  interrupt  request  on  the  rising  or 
falling  edge  of  the  input  applied  by  the  peripheral  device.  The  other  two 
control  lines  may  be  programmed  to  act  as  either  input  or  output.  In  the 
present  program  configuration,  they  are  both  used  as  outputs  to  strobe  the 
two  translation  stage  drivers.  Presently,  the  interrupt  request  feature  of 
the  PIA  is  not  utilized.  The  microcomputer  operates  at  a basic  cycle  time 
of  2 microseconds. 

Microprocessor  Controlled  Program 

The  goal  was  to  develop  a program  which  would  control  the  stepping  velocity 
and  number  of  steps  of  both  the  X-axis  and  Y-axis  stepping  motors  on  the 
Aerotech  translation  stages.  This  would  provide  the  control  needed  to  pro- 
vide laser  beam  exposure  such  as  the  LHSl  optical  waveguide  using  the  argon 
laser. 

The  approach  taken  in  designing  the  software  program  involved  building  a set 
of  short  subroutines  which  control  the  functions  in  the  translation  stages, 
configure  the  PIA  control  registers,  perform  double  word  arithmetic,  deter- 
mine time  delays  for  pulsing  the  axis,  count  the  number  of  pulses,  input 
initial  parameters  from  the  TTY,  convert  decimal  input  to  binary,  control 
the  laser  shutter,  and  supply  a debug  subroutine  which  will  output  the  contents 
of  the  program  counter,  X-register,  A-register  and  B-register. 

The  program  in  it's  present  form  contains  subroutines  which  are  not  utilized 
but  available  for  future  improvements  in  the  program.  The  program  has  not 
been  optimized  with  respect  to  execution  time  and  also  contains  many  "no 
operation  (NOP)"  commands  which  serve  as  break  points  for  insertion  of  a 
"jump  to  debug"  subroutine.  The  NOP  command  also  allow  temporary 
modification  of  the  program  without  the  necessity  of  generating  a new  object 
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tape  on  the  Honeywell  computer  network  (IlCN)  which  is  a very  time  consuming 
process.  As  the  number  of  modifications  accumulate,  it  eventually  becomes  | 

necessary  to  return  to  the  HCN  system  to  update  the  source  file  program  and 
generate  a new  assembly  listing  and  object  tape.  . 

Figure  21  shows  the  main  program  flow  chart  and  Figure  22  shows  the  Pulse 
Command  word  modification  subroutine  flow  chart.  The  present  program 
begins  by  allowing  entry  of  initial  decimal  parameters.  The  parameters  are  ^ 

converted  to  binary  form  and  scored  in  constant  locations.  They  consist  of  a * 

step  count  along  the  X-axis  (used  to  write  horizontal  line),  a step  count  which 
is  used  in  the  curve  subroutine  (counts  X-axis  pulses  only)  and  a step  count 
which  is  used  in  the  diagonal  subroutine  (again  only  X-axis  pulses  are  counted). 

Multiplier  constants  are  entered  for  use  in  the  curve  subroutine.  These 
multipliers  determine  the  number  of  pulses  issued  at  any  particular  rate 
before  the  rates  are  changed  by  fixed  delta  constants  previously  entered  for 
the  X and  Y axes.  The  multipliers  essentially  determine  the  radius  of  the 
curve.  The  X and  Y axes  stepping  rate  will  also  be  a factor  in  radius  deter- 
mination. 

When  the  X count  for  the  horizontal  line,  curves  line,  and  diagonal  have  been 
entered  and  converted  to  binary  form,  the  program  allows  entry  of  the  X and 
Y axes  multipliers,  the  initial  X and  axes  velocities  in  the  form  of  words 
to  be  counted  down  and  the  fixed  delta  counts  which  will  be  added  to  the  X 
count  thus  decreasing  the  X velocity  by  the  fixed  delta,  and  subtracted  from 
the  Y count  thus  increasing  the  Y velocity  by  the  fixed  delta.  These  decimal 
entries  are  also  converted  to  binary. 

After  the  parameters  are  entered,  the  program  proceeds  to  configure  the 
PIA,  One  eight  bit  port  interfaces  the  X-axis  driver  and  the  second  eight- 
bit  port  interfaces  the  Y-axis  driver.  The  data  direction  mask  for  each 
port  is  written  in  the  appropriate  data  direction  register  in  the  PIA  then  the 
I/O  capability  becomes  active  upon  the  command  to  open  the  peripheral  data 
r registers. 
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Figure  21.  LBES  Control  Main  Program  flow  chart.  (1  of  5) 
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P’igure  21,  [,BKS  Control  Main  Program  flow  chart.  (2  of  5) 
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Figure  21.  LBES  Control  Main  Program  flow  chart.  (3  of  5) 
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Figure  22.  I.RES  control  pulse  command  word  modification  flow  chart. 


55 


Next  the  port  words  are  initialized.  When  initialized  the  following  conditions 
are  set  up.  All  command  bits  are  set  to  zero,  then  bit  1 is  set  high.  This 
procedure  resets  the  "home  select"  F/F  in  the  driver.  On  the  X-axis  port, 
bit  5 is  held  high  to  keep  the  laser  shutter  in  the  closed  position. 

When  the  ports  have  been  initialized,  the  program  will  now  loop  until  the  home 
select  switch  is  actuated.  When  the  switch  sets  the  "home  select"  F/F,  bit  0 
of  the  status  word  goes  high  and  causes  the  program  to  exit  the  wait  loop  and 
enter  the  main  program. 

In  the  main  program  the  first  step  is  to  set  up  the  X and  Y axes  drive  direction. 
In  this  case  X will  drive  clockwise  and  Y will  drive  counter-clockwise.  The 
laser  shutter  is  then  opened,  the  horizontal  count  is  set  up  in  the  counter  and 
the  initial  pulse  command  word  is  entered.  The  horizontal  pulse  loop  is  now 
entered  and  will  not  exit  until  the  total  X count  equals  the  required  X count. 
When  in  the  horizontal  pulse  loop,  the  pulse  command  word  is  not  modified 
by  the  delta  X count. 

Upon  completion  of  the  horizontal  count,  the  program  jumps  to  the  curve 
subroutine.  A new  count  is  entered  in  the  X counter.  Both  AX  count  and 
AY  count  are  allowed  to  modify  the  pulse  command  word  when  the  multipliers 
have  counted  to  zero.  AX  count  is  added  to  the  X pulse  command  word  which 
causes  the  pulse  command  word  to  increase  in  size  thereby  requiring  more 
time  to  count  to  zero  and,  consequently,  the  velocity  along  the  X-axis,  V^, 
becomes  slower.  AY  count  is  subtracted  from  the  Y pulse  command  word 
which  causes  the  word  to  decrease  in  size  thus  requiring  less  time  to  count 
to  zero  and,  consequently,  the  velocity  along  the  Y-axis,  Vy,  becomes 
faster.  As  an  initial  experiment,  we  have  fixed  both  AX  count  and  AY  count. 

When  the  X counter  reaches  the  required  X count,  the  program  branches  to 

the  diagonal  subroutine  which  essentially  freezes  the  X and  Y pulse  command 

words  thus  fixing  the  and  . The  resulting  diagonal  continues  until  the 

A Y 


diagnonal  X count  reaches  the  required  X count  entered  in  the  initial  para- 
meters. The  program  exits  the  diagonal  loop,  the  laser  shutter  is  closed, 
and  the  computer  is  returned  to  MIKHUG  control, 

A number  of  channel  waveguides  based  on  the  laser  heating  strain  induced 
(LHSl)  technique  was  made  using  the  LBES  facility.  The  results  are  dis- 
cussed in  the  following  subsection. 


PROPERTIES  OF  WAVEGUIDE  CIRCUIT  ELEMENTS 

Demonstration  of  some  of  the  simple  waveguide  circuits  that  can  be  made 
using  the  LHSI  technique  were  performed.  In  the  following  we  discuss  the 
results  of  the  properties  of  these  waveguide  devices. 


LHSI  Waveguide  Constriction 

A continuous  straight  LHSI  waveguide  consisting  of  a wide  part  and  a narrow 
part  joined  together  was  made  by  power  reduction  to  8 0%  of  the  initial  value 
in  the  middle  of  a scanning.  The  resultant  waveguide  has  a cross-sectional 
dimension  of  lOum  radius  in  the  wide  part  and  a 15um  radius  in  the  narrow 
part.  The  waveguide  acceptance  angle  were  then  measured  for  both  the 
narrow  and  the  wide  part  of  the  constriction  with  beam  input  coupled  from 
the  wide  waveguide.  The  result  in  shown  in  Figure  23.  It  is  most  unexpected 
to  observe  that  the  acceptance  angle  for  the  TM  and  TE  modes  for  both  parts 
of  the  constricted  waveguide  to  be  about  the  same  minutes).  This  result 
may  be  due  to  the  efficient  coupling  from  the  wide  guide  to  the  narrow  guide 
by  the  tapered  constriction. 


PRISM  INPUT  ANGLE 


Figure  23.  Relative  beam  intensity  measured  at  narrow  and  wide  part 
of  a constructed  LHSI  waveguide,  as  a function  of  prism 
input  angle.  The  data  for  the  wide  part  are  • - TM  and  0 - 
TE  mode  and  those  for  the  narrow  part  are  X - TM  and 
A - TE  modes. 
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LHSI  Waveguide  Bends 


The  waveguide  bending  loss  could  be  treated  by  an  analysis  given  by  Gloge  . 
To  be  useful  in  constructing  integrated  optical  circuits,  LIISl  guides  will  have 
to  possess  acceptably  low  losses  in  the  bends  used  to  distribute  the  light 
along  the  various  branches.  We  consider  these  losses  here,  again  using  the 
parabolic-index  fiber  as  our  model. 

The  major  losses  in  multi-mode  fibers  are  caused  by  higher-order  modes 
striking  the  interface  at  less  than  the  critical  angle.  These  modes  become 
leaky  rather  than  bounded  and  after  some  distance  have  lost  their  energy. 

The  fraction  of  modes  lost  in  a bend  of  radius  R in  a parabolic-index  fiber 
. 15 

IS 

„ - 2a  (39) 

^ RA 

where  a is  the  fiber  radius,  A = ^ ^b^  ~ 2An/n^  and  An  = n^  - n^. 

This  is  independent  of  the  angle  of  the  bend,  and  is  actually  an  upper  limit 
on  loss  assuming  all  energy  in  a leaky  mode  is  immediately  lost. 

Typical  LHSI  guides  have  a = 50um,  An  = lO"^,  and  n^^  = 1.  5,  which  gives 
A = 1.  33  X lO”^  and  the  product  yR  = 7.  5cm.  So  to  limit  the  fraction  of  mode 
becoming  leaky  to  y = 0.  1 (0.  5dB)  requires  a bend  of  75cm. 

Those  bending  losses  and  minimum  bend  radii  are  probably  acceptable  for 
the  construction  of  integrated  optical  circuits,  although  it  may  be  desirable 
to  obtain  guides  with  tighter  confinement  (larger  An/n)  to  lower  these  losses. 

Experimentally,  we  have  made  waveguide  bends  by  the  LHES,  with  a bending 
radius  approaching  10cm.  Figure  24(a)  shows  such  a waveguide  bend,  fig- 
ure 24(b)  depicts  a HeNe  laser  beam  proagating  though  such  a bend.  The 
beam  intensity  measured  along  the  bend  gives  the  bending  loss  of  5db/cm. 
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Figure  24.  Photograph  of  a LliSl  waveguide  bend  produced  by  LBES; 

(a)  dark-field  photograph  of  waveguide;  (b)  photograph 
of  laser  beam  propagating  through  the  bend. 
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LHSI  Waveguide  Gaps 


We  have  also  made  straight  waveguides  with  a gap  of  various  length.  The 
gap  should  serve  as  a mode  selector  since  higher  order  modes  leaving  the 
end  of  a straight  waveguide  will  not  be  able  to  couple  into  another  section  of 
the  waveguide  at  a distance  away.  This  will  result  in  a reduction  of  the  beam 
intensity  in  the  ITISI  waveguide  as  the  gap  width  increases.  We  have  observed 
this  effect  by  measuring  the  intensity  of  the  beam  at  the  output  prism  coupler 
with  the  input  prism  at  a fixed  distance  (2.  3cm)  away  from  it.  This  2.  3cm 
waveguide  path  includes  a gap  of  varying  width.  An  example  of  the  measured 
output  is  shown  in  the  following  tabulation  (arbitrary  units  of  intensity): 


Distance 
Measured  from 
Input  Prism 


1.  5cm 

2.  3cm 


Gap  Spacing  Beginning  2cm  from  Input  Prism 
0 0,  5mm  2,  5mm 

5,42  3.  92  3.  29 

3.  13  1.  52  0.  163 


It  is  seen  that  the  additional  loss  is  1.  7db  for  the  0.  5mm  gap  and  lldb  for  the 
2.  5mm  gap.  Giving  an  average  loss  of  4db/mm  of  gap  spacing. 


LHSI  Waveguide  Coupler 

A waveguide  directional  coupler  was  designed  based  on  the  fact  that  optical 
wave  can  propagate  through  a gap  and  a bend  in  the  LHSI  waveguide  as  dis- 
cussed previously,  A schematic  diagram  of  such  a coupler  is  shown  in  Figure 
25,  The  wave  propagating  from  port  A will  exist  at  H through  a bend  and  at 
C through  a gap.  But  there  is  very  little  coupling  to  port  D if  coupling  between 
adjacent  waveguides  is  small.  Similar  result  is  obtained  for  input  ^rom  C. 


Figure  25.  Schematic  of  a I.HSI  waveguide  directional  coupler. 


To  assure  that  coupling  between  adjacent  waveguides  are  small,  we  have 
made  an  analysis  of  the  coupling  between  modes  of  adjacent  multimode  para- 
bolic-index fibers  which  can  be  applied  to  IJISI  guides.  The  mathematics 

1 6 

are  involved  and  follow  the  method  of  Marcuse^  . We  consider  the  modes  of 
each  guide  to  be  unaffected  by  the  presence  of  the  other  with  the  electric  field 
amplitudes  given  by  the  same  expressions  used  in  the  preceeding  section. 

Some  integrals  must  be  evaluated  numerically. 

The  results  of  those  calculations  are  best  expressed  qualitatively  as  follows: 

• For  a mode  described  by  mode  indices  m,  n in  Cartesian  coordi- 
nates, coupling  is  limited  to  modes  in  the  other  guide  with  indices  m', 

I 

n'  satisfying: 

it 

r, 

•> 
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m',  n' 


m,  n 

m+2,  n-2 
m-2,  n+2 


(40) 


Since  mode  effective  index  depends  only  on  the  sum  (mfn+1),  we 
see  that  the  modes  which  are  coupled  propagate  with  the  same 
phase  velocity. 

• For  a fixed  spacing  between  adjacent  guides,  only  the  highest- 
order  modes  (lowest  effective  index)  which  can  still  be  identified 
as  belonging  to  a single  fiber  are  significantly  coupled.  Lower- 
order  modes  are  much  more  tightly  confined  by  the  graded  index 
profile,  overlap  very  little,  and  are  coupled  more  weakly  by 
several  orders  of  magnitude. 

• For  a given  pair  of  coupled  modes,  the  strength  of  the  coupling 
decreases  very  quickly  with  increasing  guide  separation,  typically 
several  orders  of  magnitude  for  a change  in  separation  of  one 
diameter. 

In  summary,  appreciable  coupling  occurs  only  between  a few  modes  at  a 
time  in  adjacent  multimode  parabolic-index  guides.  This  renders  unsuitable 
the  couoled-adjacent-waveguide  approach  used  with  single-mode  step-index 
guides  or  designing  directional  couplers,  etc.  We  are  therefore  assured 
that  as  long  as  the  parallel  section  of  the  directional  coupler  is  more  than 
one  radius  of  the  waveguide  diameter  apart,  cross  talk  will  be  negligible. 

We  have  made  one  such  directional  coupler  with  a gap  scpaing  of  1.0mm  and 
a bending  of  5 degrees.  The  directionality  of  the  coupler  was  easily  observed. 
However,  the  losses  of  the  coupler  is  about  lOdb  to  each  forward  branch. 
Further  improvement  is  needed  to  reduce  this  loss. 
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SECTION  IV 


RESIDUAL  AND  INDUCED  STRAIN  IN 
WAVEGUIDE  MATERIALS 

POLISHING-INDUCED  RESIDUAL  STRAIN 

When  a surface  is  abrasively  machined  by  finish  grinding,  lapping  and  polish- 

ing,  a highly  deformed  layer  adjacent  to  the  surface  is  produced  . In  non- 

metallic  crystalline  solids  with  high  yield  points  the  deformed  layer  is  fairly 

shallow  and  the  boundary  between  the  deformed  zone  and  the  underlying  bulk 

17-20 

material  can  be  very  sharp  . The  residual  stresses  (and  strains)  in 

these  mai  hined  layers  are  compressive  if  extensive  heating  does  not  occur 
20 

during  machining  . In  polishing  operations  heating  is  minimal.  The  resi- 
dual compressive  stresses  arise  from  the  fact  that  the  deformation  taking 
place  during  polishing  is,  in  essence,  a type  of  "smearing"  that  tends  to 
increase  the  lateral  dimensions  of  the  polished  surface.  The  increase  in  sur- 
face area  is  resisted  by  the  underlying  bulk  and  the  deformed  layer  is  put  in 
a state  of  biaxial  compression.  If  the  strain  optic  coefficient  where  n is 
the  index  of  refraction  and  S is  the  strain,  is  negative,  the  index  in  the 
deformed  layer  will  increase  as  a result  of  polishing  and  the  layer  could  act 
as  a waveguide.  This  is  the  situation  found  by  Chen  et  al  in  polished  PLZT^, 

In  this  study  we  have  extended  their  observations  on  PLZT  and  have  syste- 
matically attempted  to  introduce  a waveguiding  surface  layer  in  a material 
with  a negative  strain  optic  coefficient  and  a much  simpler  microstructure 
than  PLZT,  namely  single  crystal  lithium  niobate  (LiNbO^).  In  what  fol- 
lows, the  results  of  the  study  of  PLZT  waveguides  are  presented  first 
followed  by  the  discussion  of  waveguiding  in  LiNbO^  (LN)  crystals. 
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Waveguiding  in  PLZT 

In  a material  such  as  PLZT  ceramic,  the  magnitude  of  the  photoelastic  coef- 

25 

ficient  p is  expected  to  be  about  0.  2 . The  polishing- induced  strain  S for 

17  - 3 - 2 

this  case  could  be  in  the  range  of  10  <S<10  . Since  the  refractive 

_ 2 

index  is  2.  5,  if  we  assume  S = 10  , the  polishing- induced  refractive  index 

”2  1 8 
could  be  as  high  as  10  . If  the  layer  thickness  W is  a few  pm,  as  expected 

there  will  be  a few  TE  and  TM  modes  that  could  propagate  in  the  polished  layer 

of  PLZT  wafer.  Experimentally,  we  have  achieved  optical  waveguiding  in 

polished  PLZT  wafer.  The  measured  waveguide  effective  refractive  index 

for  the  two  lowest-order  TE  modes  gives  a best  fit  value  of 

An  ~ 0.  018  and  W ~ 2.  33pm 

This  result  was  obtained  in  an  earlier  research  work  entitled  "An  Investiga- 
tion of  PLZT  for  Interpreted  Optics"  supported  by  the  Office  of  Naval 
Research.  Further  details  could  be  found  in  the  final  report^  of  that  research 
work. 


Polishing- Induced  Strain  in  LN  Crystals 

Our  choice  of  studying  strain-induced  waveguiding  in  a single  crystal  as 
opposed  to  a polycrystalline  material  such  as  PLZT  is  based  on  the  idea  that, 
first,  characterization  of  the  waveguiding  layer  is  easier  in  the  crystal  and, 
second,  that  the  anisotropic  properties  of  the  crystal  could  be  used  to  our 
advantage.  Briefly,  the  study  consisted  of  first  giving  different  crystal- 
lographic faces  a number  of  different  polishing  treatments.  The  nature  and 
extent  of  deformation  and  damage  in  the  polished  layers  were  then  charac- 
terized by  x-ray  techniques  and  by  direct  observation  in  an  optical  and  a 
s canning  electron  microscope.  Finally,  the  polished  surfaces  were 
then  checked  to  determine  if  they  guided  a wave  from  a HeNe  laser. 
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Experimental  Procedure  --  The  experimental  procedure  employed  is  des- 
cribed in  the  following  paragraphs. 

Polishing  Technique  --  The  samples  were  polished  in  a five-step  opera- 
tion similar  to  that  used  to  prepare  surfaces  for  microscopic  examination. 
For  lapping  and  polishing  the  samples  were  mounted  on  metal  discs  with  a hot 
melting  wax.  Cylindrical  conditioning  rings  were  used  to  maintain  flatness 
during  lapping  and  lithium  niobate  "feet"  were  cemented  to  the  outside  of  the 
specimen  holder  to  maintain  flatness  during  polishing.  The  five  stages  of 
surface  preparation  were; 

• Stage  1 --  Lapped  with  a 12pm  SiC  slurry  on  cast  iron  using  a Gyro- 
matic  12  precision  lapping  machine.  An  abrasive  concentration  of 
20gm /liter  of  distilled  water  was  used  at  a rotation  rate  of  60  rpm. 
When  the  downward  stress  on  the  sample  was  3 psi  the  material 
removal  rate  was  approximately  12pm /min. 

• Stage  2 --  Lapped  with  a 3pm  AI2O3  slurry  on  cast  iron  using  a 
Frantz  "unipol"  lapping  machine.  An  abrasive  concentration  of 
lOgm/liter  of  distilled  water  was  used  at  a rotation  rate  of  90  rpm. 
The  removal  rate  was  approximately  6pm/min. 

• Stage  3 --  Lapped  for  15  minutes  with  the  same  AI2O3  slurry  on 
the  same  machine  using  a lead  lap.  The  removal  rate  was  ap- 
proximately 1pm /min, 

• Stage  4 --  Hand  polished  for  several  minutes  with  3pm  diamond 
abrasive  on  Metcloth*  rotating  at  550  rpm.  Kerosene  was  used 
as  a carrier  for  the  abrasive. 


*A.  B.  Buhler  Co. 
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• stage  5 --  Polished  for  one  to  two  hours  with  cerium  oxide  on  a 

Pellon  PANW  lap  using  a Strasbough  RGY  Precision  Polishmaster 
machine  operating  at  8 0 rpm. 

Examination  of  Lapped  and  Polished  Surfaces  --  The  LiNbO^  samples 
were  examined  at  all  stages  of  preparation  optically  and  with  a Cambridge 
Instruments  Model  S4-10  scanning  electron  microscope.  To  prevent  charging 
while  in  the  SEM  the  samples  were  coated  with  carbon. 

Characterization  of  Surfaces  by  X-Ray  Techniques  --  One  of  the  ways  to 
characterize  the  highly  deformed  subsurface  layer  introduced  by  lapping  and 
polishing  is  by  x-ray  techniques.  The  residual  stresses  in  the  layer  cause 
changes  in  the  lattice  parameters  which  may  be  detected  by  x-rays.  The 
changes  may  be  subtle  so  precision  techniques  are  necessary.  In  this  study 
we  used  Bond's  technique  to  measure  the  changes  in  lattice  parameter  per- 
pendicular to  the  surface  resulting  from  the  biaxial  stresses  in  the  layer 
parallel  to  the  surface.  The  two  are  related  by  Poisson's  ratio  with  the 
changes  perpendicular  to  the  stresses  being  about  three  times  the  changes 
parallel  to  the  surface. 

The  Bond  technique  uses  a symmetrical  diffractometer  sue'  as  shown  in 
Figure  26.  The  crystal  is  rotated  till  a reflection  is  detec-ad  in  the  upper 
detector.  It  is  then  rotated  till  the  same  reflection  is  detected  in  the  lower 
detector.  By  accurately  measuring  the  angle  between  the  two  reflecting 
positions  the  lattice  parameters  can  be  measured  to  a few  ppm. 

The  necessary  sensitivity  can  be  discussed  in  terms  of  the  angular  changes 
that  need  to  be  measured  to  detect  strains  in  the  deformed  layer  on  the  order 
of  1 to  10  ppm.  The  sensitivity  requirements  for  x-ray  -asurements  of 
lattice  parameters  are  found  from  Bragg's  Law,  nX  = 2d  sin  fi  where  n is  an 
integer,  X the  wavelength,  d the  interplanar  spacing  and  0 the  diffraction 
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angle.  Straining  the  lattice  causes  a change  in  d.  Thus,  from  the  first 
derivative  of  Bragg's  Law 


A9  = - ^ tan  e (41) 

d 

we  see  that,  for  a constant  strain  (S  = the  corresponding  change  in 

diffraction  angle,  A 9,  is  proportional  to  the  tangent  of  the  diffraction  angle. 

Since  tan  9 increases  rapidly  near  9 = 90°,  most  high-sensitivity  lattice 

parameter  measurements  are  made  in  this  region.  The  increase  in  A 9 as 

the  diffraction  angle  approaches  90°  is  shown  in  Figure  27  for  lattice  strains 
- 1 -6 

of  10  and  10  . Values  of  AOas  small  as  5 to  10  arc  seconds  can  be 

detected  by  our  Bond  diffractometer.  As  discussed  below  diffraction  angles 

of  around  75  degrees  were  used  to  investigate  LN.  As  seen  from  Figure  27 

- 5 

this  enabled  us  to  detect  lattice  strains  on  the  order  of  10  ‘ with  the  dif- 
fractometer. A brief  description  of  the  diffractometer  is  given  next  followed 
by  a discussion  of  the  accuracy  and  errors  inherent  in  the  device. 

The  Bond  symmetrical  diffractometer,  as  shown  in  Figure  26,  consists  of 
three  main  components: 

• An  x-ray  source  and  collimator 

• A specimen  holder  attached  to  a sensitive  angle  measuring  device 

• Dual  x-ray  detectors 

The  diffraction  angle  9 is  obtained  from  9°(deg)  = 90°  - (N^-N2)/2  where 
and  are  the  angular  positions  of  the  diffracting  plane  normals  in  the  upper 
and  lower  diffracting  positions  respectively.  Since  the  highest  resolution  is 
obtained  at  diffraction  angles  near  90  degrees  the  x-ray  source  must  be  prop- 
erly selected. 
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UPPER 

detector 


D(FFRA(^TIOI\l  ANGLE  9 = 90  - ^/z 
N,  ANDN2ARE  DIFFRACTION  PLANE 

normals  at  upper  and  lower 

DIFFRACTING  POSITIONS  RESPECTIVELY 


Figure  27.  Change  in  diffraction  angle  A0  caused  by  strain 
plotted  as  a function  of  diffraction  angle,  B. 
Curves  for  strains  of  lO"^  and  10~6  are  plotted. 
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The  symmetrical  configuration  eliminates  many  errors.  However,  some 
errors  such  as  incident  beam  axial  divergence,  refraction  and  the  Lorentz 
Polarization  Factor  error  are  still  present.  These  errors  can  be  readily 
treated  mathematically  and,  in  fact,  are  small  at  large  values  of  0. 

The  x-rays  are  produced  by  a standard  x-ray  tube  with  a 1mm  x 10mm  line 
source.  Slits  positioned  between  the  source  and  the  sample  allow  a 150pm 
X 1mm  beam  to  impinge  on  the  sample.  The  crystals  were  mounted  on  the 

22 

rotating  sample  holder  of  a Phillips  powder  diffractometer.  Beu  has  shown 

the  angular  sensitivity  of  this  type  of  support  is  insufficient  to  measure  the 

- 5 

A0of  5 to  10  seconds  needed  to  detect  10  strain  (see  Figure  27).  For  this 

reason  an  autocolimator  and  high  quality  optical  angle  standards  were  used  to 

make  the  measurements.  This  is  a standard  technique  for  making  angle  mea- 
23 

surements  . With  this  setup  we  were  able  to  measure  A0to  4 seconds. 

Results  and  Discussion  --  Experimental  results  are  discussed  in  the  follow- 
ing paragraphs. 

SEM  Characterization  of  Lapped  and  Polished  LN  Crystal  Faces  --  Our 
goal  in  this  portion  of  the  study  was  to  introduce  plastically  deformed  layers 

24 

on  LN  crystal  faces  by  lapping  and  polishing.  Based  on  previous  experience 
the  plastic  deformation  resulting  from  an  abrasive  machining  operation  can 
be  revealed  by  electron  microscopy.  In  the  present  work,  scanning  electron 
micrographs  were  made  of  the  LN  crystal  faces  after  each  lapping  and 
polishing  step  to  follow  the  development  of  the  deformed  layer.  Figures  28 
through  31  are  scanning  electron  micrographs  at  2000X  which  illustrate  the 
nature  of  tlie  four  crystal  faces  after  each  of  the  five  lapping  and  polishing 
steps.  Figure  28  shows  the  effects  of  machining  on  the  (0001)  face.  The 
rough  regions  on  the  lapped  surfaces  are  where  material  has  been  removed 
by  brittle  fracture.  The  smooth  patches  are  where  material  was  removed 
by  plastic  deformation,  i.  e.  , the  regions  are  "smeared".  As  finer  lapping 
conditions  are  used,  the  extent  of  brittle  fracture  on  the  surfaces  decreases. 
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3u  Diamond 


Figure  28.  Scanning  electron  micrographs  showing  surface  structure  of 
(0001)  LN  face  following  successive  lapping  and  polishing 
operations. 


3u  Diamond 


Figure  29.  Scanning  electron  mi''rographs  showing  surface  structure 
of  (1010)  ].N  face  following  successive  lapping  and 
polishing  operations. 


Figure  30.  Scanning  electron  rnicrographs  showing  surface  structure  of 
(lOil)  LN  face  following  successive  lapping  and  polishing 
operations. 
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Figure  31.  Scanning  electron  micrographs  showing  surface  structure  of 
(1120)  l.N  face  following  successive  lapping  and  polishing 
operations. 
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After  polishing  with  3pm  diamond  all  indications  of  fractured  regions  are 
gone,  the  surface  is  smooth  and  contains  fine  scratches  left  by  the  abrasive. 
After  a final  CeO  polish  on  Pelion  the  scratches  are  gone  and  the  surface 
is  smooth.  The  surface  features  on  the  micrograph  of  the  polished  surface 
presumably  are  polishing  debris  and  imbedded  abrasive  particles.  Material 
removal  during  the  production  of  the  last  two  surfaces  was  clearly  by  plastic 
deformation.  The  extent  of  the  plastically  deformed  zone  and  the  residual 
stress  state  will  be  discussed  later. 

The  stepwise  development  of  the  polished  surface  on  a (1010)  LN  crystal 
face  is  shown  in  Figure  29.  The  increasing  amount  of  plastic  deformation  as 
the  lapping  treatment  becomes  finer  is  again  obvious.  No  differences  between 
the  (0001)  and  (lOTO)  surface  polished  with  3pm  diamond  and  with  CeO 
are  evident.  The  size  of  the  fractured  zones  on  the  lapped  (1010)  faces 
appear  slightly  larger  than  those  on  the  lapped  (0001)  faces.  This  could,  in 
part,  be  due  to  the  fact  that  the  (10T2)  cleavage  plane  makes  a 16-degree 
angle  with  (lOTO)  while  the  angle  between  the  cleavage  plane  and  the  basal 
plane  is  about  60  degrees. 

A similar  series  of  scanning  electron  micrographs  after  lapping  and  polish- 
ing operations  have  been  performed  on  (lOTl)  faces  is  shown  in  Figure  30. 

The  3pm  diamond  polished  surface  showed  features  that  appear  to  be  pits. 

The  last  micrograph  in  the  series  show,  the  depressions  were  not  totally 
removed  by  final  polishing  in  this  case.  The  (10T2)  cleavage  plane  makes 
an  angle  of  16  degrees  with  the  (1011)  face  which  may  again  account  for  the 
larger  patches  on  the  lapped  surfaces  where  material  was  removed  by  brittle 
fracture. 

Figure  31  shows  the  (11^0)  crystal  face  after  various  treatments.  The 
straight,  flat  regions  in  the  first  micrograph  show  where  material  was 
"smeared"  by  the  large  abrasive  particles.  The  regions  where  material 
has  been  removed  by  brittle  fracture  appear  smaller  in  this  case  as  expected 
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since  the  cleavage  plane  makes  an  angle  of  55  degrees  with  the  machined  sur- 
face. Some  pits  are  present  after  3pm  diamond  polishing  which  persist  after 
polishing  with  CeO.  The  surfaces  with  pits  and  scratches  are  clearly  not 
acceptable  to  use  as  waveguides  because  the  imperfections  would  scatter  the 
beam.  The  polishing  techniques  used  on  these  surfaces  would  have  to  be 
refined.  It  is  interesting  that  the  residual  damage  in  the  (1011)  and  (1120) 
polished  surfaces  is  higher  than  in  the  (0001)  and  (1010)  surfaces  since  all 
were  polished  simultaneously. 

Precision  Lattice  Parameter  Measurements  --  The  lattice  parameters 
measured  on  the  LN  crystal  faces  lapped  and  polished  in  this  study  are  listed 
in  Table  3.  As  noted  from  the  table,  measurements  were  made  after  each 
surface  finishing  operation.  Also  included  on  the  table  are  lattice  parameter 
measurements  made  on  commercially  polished  (0001)  and  (1010)  faces  and  on  a 
(lOTO)  face  with  a Ti  in-diffused  layer  on  it.  The  precisions  indicated  on  the 
table  are  estimates  based  on  a measurement  capability  of  ±4  arc  seconds 
with  our  autocollimator  and  the  following  diffraction  angles: 


Diffraction 


Crystal  F’ace 

Plane 

0 (deg) 

Tube 

(0001) 

(00018) 

74.  5 

(lOTO) 

(0060) 

8 2.  5 

(llIO) 

(3360) 

70.  4 

CoKg 

(lOTl) 

(6066) 

8 0.  1 

CuKg 

In  every  case  the  surface  treatment  caused  an  increase  in  the  lattice  spacing. 
It  is  the  Poisson  expansion  due  to  the  compressive  stresses  in  the  deformed 
*■  surface  layer  introduced  by  the  machining  treatment.  A better  indication  of 

the  magnitudes  of  the  changes  in  lattice  spacing  is  given  in  Table  3 which 
lists  the  percentage  change  in  spacing,  i.  e.  , the  strain.  The  lowest  strain 
was  observed  on  the  Z-cut  (0001)  face.  The  strain  on  the  x-cut  (lOTO)  face 
did  not  change  systematically  with  treatment.  The  strain  on  the  (1010)  face 
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Table  3.  Effect  of  Surface  Treatment  on  Interplanar 
Spacing*  of  Lithium  Niobate  Crystals 


Machined  Crystal  Face 

Surface  Treatment 

c 

(0001) 
±.  0001 

a 

(1010) 

±. 00002 

b 

(1120) 

±.  00002 

(1011) 

±.  00005 

Annealed** 

I2u  SiC  on  cast  iron  lap 
3u  Al^O,^  on  cast  iron  lap 
3u  Al^O^  on  lead  lap 
3u  diamond  on  Metcloth 
CeO  on  Fellon 


Commercially  polished 


Ti  diffused 


13.863 
13.  8647 
13. 8645 
13. 8647 
13.  8647 
13. 8649 


5.  1483 

5.  1483 

5. 1498 

5.  1501 

5. 1499 

5.  1508 

5. 1499 

5.  1501 

5. 1500 

5.  1491 

5.  1501 

5.  1499 

5. 1500 

5.  1499 

5.  1499 

4. 2444 
4. 2459 
4. 2456 

4. 2458 

4. 2459 
4. 2458 


* These  values  are  not  corrected  for  refraction  Lorentz  polarization  factor 
or  for  axial  divergence  of  the  incident  beam. 

'■*W.  L.  Bond  referenced  by  S.  C.  Abrahams,  et.  al. , J.  Phys,  Chem.  Solids, 
27,  997  ( 1966). 
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^ Table  4.  Percentage  Change  in  Lattice  Parameter  (i.  e. , strain) 

j as  a Function  of  Surface  Treatment  of  LN  Crystals 


Surface  Treatment 

Z Cut 
(GOOD 

X Cut 
(1010) 

Y Cut 
(1120) 

— 

(loll) 

12w  SiC  on  cast  iron  lap 

.0123 

o 

to 

CD 

. 0350 

. 0353 

3u  A1„0„  on  cast  iron  lap 

. 0108 

.0311 

. 0486 

. 0283 

3u  A1„0^  on  Pb  lap 

. 0108 

. 0311 

. 0350 

. 0330 

3iu  diamond  on  Metcloth 

. 0123 

. 0330 

. 0155 

. 0353 

CeO  on  Pellon 

. 0137 

. 0350 

.0311 

. 0330 

Commercially  polished 

. 0123 

. 0330 

. 0311 

Ti  diffused 

.0311 

I 

I 
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was  somewhat  higher  and  did  not  appear  to  increase  as  the  surface  treatment 
produced  a finer  surface.  Results  obtained  on  the  Y-cut  (1120)  and  (lOTl) 
faces  showed  strains  about  as  large  as  those  observed  on  the  (lOTO)  face.  No 
systematic  change  with  treatment  was  observed.  Results  on  the  commercially 
polished  samples  showed  the  strains  to  be  the  same  as  those  produced  by  us 
with  3^m  diamond  or  with  CeO.  The  strain  in  the  Ti  diffused  (lOTO)  face  is 
slightly  less  than  that  produced  by  polishing. 

Optical  Waveguiding  in  Polished  Samples  --  After  fabrication  the  polished 
crystal  surfaces  were  evaluated  as  optical  waveguides. 

The  crystal  was  mounted  on  a rotating  head  fixture.  A 60-degree  rutile  prism 
was  used  to  couple  the  beam  into  the  polished  surface.  For  beam  polarized 
perpendicular  to  the  LN  wafer,  (TM  case)  the  coupling  angles  are  around 
-4  and  +7  degrees  for  Z-cut  and  Y-cut  wafers  respectively.  For  beam 
polarized  parallel  to  the  wafer  plane,  the  coupling  angles  are  around  -30  and  -22 
degrees  for  coupling  into  n^  and  n^  respectively.  The  crystals  are  also 
rotated  90  degrees  in  the  plane  so  that  the  beam  propagation  directions  were 
changed. 

To  the  best  of  our  ability  we  did  not  detect  any  waveguiding  in  any  of  the  sur- 
faces polished  in  our  lab  or  in  the  commercially  polished  samples.  We  con- 
clude from  these  results  that  the  residual  stresses  introduced  by  the  lapping 
and  polishing  operations  were  insufficient  to  increase  the  index  of  refraction 
in  the  deformed  layer  adjacent  to  the  polished  surface.  This  will  be  discussed 
further  in  the  following. 

Discussion  of  Results  --  Our  results  to  date  have  shown  first  that  lap- 
ping and  polishing  introduces  a measurable  amount  of  strain  into  the  surface 
of  a lithium  niobate  crystal  and  secondly  that  the  strained  surface,  contrary 
to  our  expectations,  does  not  guide  a laser  beam.  We  can  calculate  the 
change  in  the  index  of  refraction  in  the  deformed  layer  from  the  measurement 


f 
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of  strain  made  by  the  Bond  x-ray  technique.  A change  in  index  on  the  order  of 


lO"^  to  lO"^  should  be  sufficient.  The  change  of  refraction  in  the  strained 


layer  is  related  to  the  strain  by  the  fourth-rank  strain  optic-tensor 
which  couples  the  dielectric  impermeability  tensor  ( ^2  ~ 

n J J 


strain,  as  follows: 


A3-  • = P.  ■,  1 S,  , 
ij  ijkl  kl 


(42) 


Since 


B . . 

13 


The  strain-induced  change  in  index  is  then 


An.j  = -d/2)n°'^AB.^ 


(43) 


where  the  direction  of  propagation  is  perpendicular  to  X^.  and  • n°  is  the 
index  of  the  unstrained  material.  The  strain  optic  tensor  and  the  strain 
tensor  are  written  with  respect  to  the  rhombohedral  LiNbO^  axes.  The 


r LiNbO 

'3 

in  the  contracted  notation: 

Pll 

^12 

Pl3 

Pl4 

0 

0 

^21 

^22 

^23 

P24 

0 

0 

P3I 

p 

32 

^33 

0 

0 

0 

P4I 

P42 

0 

P44 

0 

0 

0 

0 

0 

0 

^55 

P| 

0 

0 

0 

0 

^65 

P, 

5R 

66 
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where 


^11 

= P22  = 

-0. 

02 

P12 

" ^21 

+0, 

08 

^33 

0 

0 

II 

P31 

^32 

+0. 

17 

Pl3 

" ^23  " 

+0. 

13 

Pl4 

= -^24  = 

= -] 

^65 

P41 

" '^42 

= P 

56 

P44 

" ^55  ^ 

+0. 

12 

■0.  08 


^66  ' ‘'2(Pn  - Pl2>  = -O' “5 


The  principal  axes  of  the  strain-tensor  lie  perpendicular  to  and  parallel  to  the 
polished  surface. 


The  values  of  the  measured  strains  perpendicular  to  the  polished  faces  are 
(from  Table  4): 


Z-cut 
X-cut 
Y - cut 


S3  = 1.  3 X 10 


-4 


= 3.  3 X 10 


-4 


83  = 3.  1 X 10 


-4 
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Consider  the  strain-induced  change  in  index  of  refraction  with  light  propa- 
gating parallel  to  the  polished  surface  of  three  LiNbO^  faces;  (1)  Z-cut  face, 
i,  e.  , (0001);  (2)  X-cut  face,  i.  e. , (1010);  and  (3)  Y-cut  face,  i.  e.  , (1120). 

(1)  Z-Cut  Face  --  In  this  case  the  strain,  S^,  perpendicular  to  the  pol- 
ished surface  is  the  component  measured  by  the  Bond  x-ray  technique; 

s = = -VS3 

V = Poisson's  ratio  = . 3 

= ^^^0  2 = P13)  S = + 1.  45  X 10  (42) 

*«3=<2P3i-ip33)S  = -. 416x10-5 
A3^  = AB5  = ABg  = 0 

(2)  X-Cut  Face  — In  this  case  the  strain,  S^,  perpendicular  to  the 
polished  surface  is  the  component  measured; 


S = $2  = S3  =vSj 

^®1  " ^ ^^11  ^Pl2^Pl3^^  " ^ 

• AB2  - <-7Pi2+Pi1+Pi3^^  = i 1.57  X 10‘ 

r AB3  = (1-(7)P3i^P33  S = +3.26x10" 

I" 

^ A64  = (-  i - DP41S  = -6.49  X 10"^ 

r 
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(3)  Y-Cut  Face  --  In  this  case  the  strain,  S^,  perpendicular  to  the 


polished  surface  is  measured: 


S = S = S S„  = -- 
1 3 2 V 


^®2  ~ ^^12"  V ^11*^13^^  ^ 

AB3  = ^^"i^P31+P33  S = -3.03x10'^ 


(44) 


A3^  = (II  = 6.  04  X lO"^ 


We  can  now  calculate  the  index  change  for  light  beams  propagating  in  the 
strained  surfaces.  There  are  12  possibilities  that  exist  depending  on  the 
propagation  direction  and  polarization.  These  are  listed  in  Table  5.  For 
simplicity  we  will  substitute  the  X^,  and  X^  axes  for  the  X,  Y,  Z axes. 


The  Z-cut  face  is  then  perpendicular  to  the  X^  axis,  the  X-cut  face  to  the 


X^  axis  and  the  Y-cut  face  to  the  X^  axis. 


The  ordinary  refractive  index  in  the  unstrained  crystal  ar  n and  n respective- 


ly. n is  in  the  X,  and  X„  directions  and  n^  is  along  the  X„  direction;  n = 2.  20, 
o Iz  0^3  o 


n = 2.  297. 
e 


(i)  face,  d = X^,  p = X^ 


An  = - 1 n ^ABo  = - 7,  68  x lO’^ 

45  O 2, 


(45) 


-> 

d = X 


1 P=^3 
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Polished  Crystal 
Face,  F 

1 

1 Propagation 

Direction,  d 

1 

Polarization 
Direction,  p 

Z-cut 

-^2 

^2 

X-cut 

^2 

^1 

^3 

^1 

^2 

V-cut 

^2 

-'^3 

'^3 

^1 

^2 

t 


■n 

f 


An=  - 1/2  n ^AP„  = 2.  52  x lO"^ 
e o 


d = P = X. 


An=  - 1/2  n ^ AB,  = - 7,  68  x lO'^ 
o 1 


d = X.,  p = X, 


An  = - l/2n  "^ABo  = < 2.  52x  lO"^ 

0 O 


(ii)  face  d = X^  p = X^ 


An-  - 1/2  n ^AB,  = + 1.  47  x 10 
o 1 


(46) 


In  this  case  the  indicatrix  ellipse  is  tilted  slightly  with  respect  to  the 
X„  and  X„  axes.  The  rotation,  however,  is  small. 


d = X2  P = X3 


3 4 

An  - - 1/2  n AP„  = - 1.  97  x lO"^ 

0 O 

d = X3  P = X^ 

~ 3 -4 

An  - -l/2n  AB,  = 1 1.47x10 
o 1 

d = X3  p = X^ 

An=  - l/2n  ^AB.,  = - 8.32x  10 
o 2 
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(iii)  face  d = p = X^ 


An--l/2n  A6.  = f 1.  36  x lO' 
o 2 


d=X^  P = X3 

An~  -l/2n  ^ABo  = - 1.81  x lO”'^ 
o 3 


d = Xg  P = 

An  = - 1/2  n ^AB.  = - 7.  42  x lO'^ 
o 1 


d = X„  p = X, 


3 -4 

An=  - l/2n^  AB^  = + 1.  36  x 10 


Note  that  in  some  cases  the  change  in  index  is  positive  and  the  magnitude  of 
~ -4 

the  change  An~10  . This  should  be  sufficient  for  waveguiding  provided  that 

the  affect  region  is  sufficiently  thick.  Calculations^  show  that  the  thickness 

of  the  layer  necessary  to  guide  a wave  depends  on  An,  Figure  32  shows  the 

lowest  order  TE  waveguide  modes  for  a substrate  of  refractive  index  n = 2.  5 

and  the  waveguide  layer  of  refractive  index  of  2.  5001,  2.  501  and  2.  51  (this 

result  is  applicable  for  n = 2.  2 and  2.  29  since  the  modification  is  negligible). 

The  coordinate  W in  nm  is  the  waveguide  layer  thickness.  In  order  to 

-4 

observe  a guided  wave,  for  An~10  , as  we  have  shown  to  be  the  case  for 

polished  LN,  the  thickness  of  the  layers  must  be  on  the  order  of  10  microns. 

The  thickness  of  the  deformed  layer  in  LN  has  not  been  measured.  However, 

the  thickness  of  the  deformed  layer  in  polished  sapphire  crystals  has  been 

1 T 

measured  by  Hockey  using  transmission  electron  microscopy  and  has  been 
found  to  be  on  the  order  of  1pm.  Presumably  the  deformed  layers  in  LN  are 
of  the  same  dimensions,  and  are  not  thick  enough  to  cause  the  strained  region 
to  act  as  a waveguide. 
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It  is  interesting  to  note  that  the  measured  lattice  strain  in  the  Ti  diffused 
layer  on  a (1010)  LN  face  is  about  the  same  as  that  of  a polished  (1010)  face. 
The  Ti  diffused  layer  acts  as  a waveguide  but  the  polished  face  doesn't.  Most 
likely  the  thickness  of  the  Ti  in  diffused  layer  is  sufficient  to  case  it  to  act  as 
a waveguide. 


MEASUREMENT  OF  PHOTOELASTIC  EFFECTS 

An  important  part  of  the  present  research  is  to  characterize  the  photoelastic 
effects  in  candidate  materials  for  the  determination  of  their  usefulness  in  the 
integrated  optics  applications.  A new  interferrometric  technique  has  been 
developed  for  the  measurement  of  this  effect.  Results  for  the  Corning  CS-2- 
60  filter  glass.  Incite,  and  LiNbOg  wafers  were  obtained.  These  as  well  as 
the  measurement  techniques  and  theory  are  given  in  the  following  paragraphs. 


Measurement  Technique  and  Theory 

From  symmetry  argument,  we  can  show  that  there  are  two  nonzero  coefficients 
in  the  photoelastic  tensor  p^^^^  and  Pj2  isotropic  material  such  as  glass. 

In  crystalline  materials,  symmetry  considerations  in  the  crystalline  structure 
will  result  in  certain  reduction  and  relationship  among  the  photoelastic 
tensor  components.  For  instance,  in  LiNbOg,  as  stated  earlier,  there  are 
only  eight  independent  nonzero  components. 

In  this  case  the  measurement  of  the  p components  will  require  crystals  of 
different  cuts.  For  our  measurements,  thin  wafers  of  crystals  of  two 
given  cuts  (X  cut  means  the  surface  of  the  wafer  is  perpendicular  to  the  X- 
axis)  were  used,  and  strain  is  applied  perpendicular  to  this  axis.  Table  6 
gives  examples  of  the  orientation  of  crystal  and  beam  for  obtaining  various 
p values  in  LiNbOg. 
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Table  6.  Crystal  and  Beam  Orientation  for  Obtaining  p Values  in  LN 


f 


[ 

i 

[ 


I 


» I 
1 


Crystal 

Cut 

p-values 

Measured 

Light  Beam 
Propagation 
Direction 

Light  Beam 

Polarization 

Direction 

Direction 
of  Strain 

Z 

Pll 

Y 

X 

X 

P3I 

Y 

Z 

X 

X 

P12 

Z 

X 

Y 

P22'"Pll 

Z 

Y 

Y 

Y 

P23=Pl3 

X 

Y 

Z 

P33 

X 

Z 

Z 

The  p's  are  indirectly  measured  by  measuring  the  refractive  index  change 

1 3 

An.  An  and  p are  related  through  An  = - n p.  S. 

To  measure  An,  a Michelson  interferometer  arrangement  was  constructed. 
The  experimental  arrangement  is  shown  in  Figure  33.  The  laser  beam  was 
first  expanded  to  a diameter  comparable  to  the  thickness  of  the  sample.  The 
final  spatial  filter  contains  a 39mm  fl  lens  to  expand  the  fringe  pattern  onto 
the  film  plane  of  a speedgraphic  camera  body. 

e 

All  measurements  were  carried  out  using  a 6328A  HeNe  laser  beam.  Let 
the  distance  for  the  reference  beam  (between  the  beam  splitter  to  mirror  M2) 
and  the  main  beam  (between  the  beam  splitter  to  Ml)  be  I'  and  f, respectively, 
and  assuming  the  physical  beam  path  length  in  the  sample  is  b,  the  phase 
difference  between  the  main  beam  and  the  reference  beam  is  given  by 

q,  = JiL  [2£>  - 2U-h)  - 2bn]  (48) 

X 
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where  n is  the  refractive  index  of  the  sample  and  X is  the  wavelength  of  the 
beam.  If  we  now  rotate  M2  around  it^  center  in  the  interferometer  plane  (x-z 
plane)  with  a small  angle  0,  then,  depending  on  the  beam  position,  we  have 

t'  = X (49) 

This  will  cause  a number  of  parallel  interference  lines  parallel  to  the  y- 
direction  to  appear  at  the  film  plane. 

To  apply  the  strain  to  the  sample  a four-point  loading  fixtured  was  used.  The 
geometry  of  the  loading  on  the  sample  is  shown  in  Figure  34. 

With  this  loading  scheme  the  stress  along  the  center  line  of  the  sample  may 

2fi 

be  calculated  from  an  elementary  beam  flexure  formula  : 


where  a = stress  along  centerline 

M = total  bending  moment  of  sample 
I = moment  of  inertia  of  sample  cross  section 
y = vertical  distance  from  neutral  axis  of  sample 


For  a sample  of  thickness  a and  width  b the  moment  of  inertia  is 


I 


ba"" 

TT 


The  total  bending  moment  exerted  on  the  sample  by  the  four  point  loading 
scheme  is  M = fc. 
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Figure  33.  Schematic  diagram  of  the  Michelson  interferrometer 
for  the  measurement  of  photoelastic  effect:  f-film, 
C“Camera,  S.  F,  -spatial  filter,  B.  S,  -beam  splitter, 
S.  H.  -sample  holder,  Ml-main  beam  mirror,  M2- 
reference  beam  mirror,  M-corner  reflector,  B.  E.  - 
beam  expander,  P-polarizer,  L-laser, 


Figure  34.  Geometry  of  four  point  loading  on  the  sample  of 
thickness  a. 


Therefore  the  stress  along  the  centerline  of  the  sample  in  the  four  point 
loading  setup  is 


a 


= 12fcy 
ba^ 


The  stress  is  directly  proportional  to  the  vertical  distance  from  the  neutral 
axis  of  the  sample. 

Therefore  the  change  in  refractive  index  along  the  centerline  of  the  sample 
is 

where  E is  the  Young's  modulus. 

With  a tilted  mirror  M2  the  interference  pattern  produced  with  an  unstressed 
sample  was  a series  of  vertical  fringes,  the  linear  refractive  index  change 
when  stress  is  applied  tilts  the  fringes  at  an  angle  to  the  vertical.  The  phase 
shift  produced  by  the  stress  induced  An  can  be  found  by  measuring  the  tilt 
angle  0 of  the  fringe  pattern  as  shown  in  Figure  35. 


Figure  35.  Tilting  of  fringe  lines  from  original  unloaded 
condition  (dotted  lines). 
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The  phase  sift  for  light  beam  pass  through  the  sample  at  distance  y from  the 
center  of  the  sample,  along  the  central  axis,  due  to  stress  is 


= 2TT  Mil 
d 

= 2bAn(y) 


(51) 


Therefore 

Comparing  expressions  for  An 

Ad  EXa^  EX  a^  , tan  0 ^ 

p.  . = ■ ■ = -I  ( — ) 

d 12fcyn3  12fcn3  d 


(52) 


(53) 


The  sample  holder  and  loading  apparatus  is  shown  in  Figure  36.  It  consisted 
of  a spring-loaded  plunger  that  could  be  screwed  down  on  a Lebow  3397  25- 
pound  load  cell  which  in  turn  exerted  pressure  on  the  upper  pressure  bars  on 
the  sample.  The  Lebow  3397  Load  Cell  is  a strain  gage  bridge  that  can  be 
biased  at  a d-c  level  to  produce  a d-c  signal  output  of  ImV  (lb-force).  The 
load  cell  was  calibrated  using  known  weights  and  a Hewlett  Packard  3476B 
digital  multimeter.  To  produce  a change  in  the  output  of  the  load  cell  of  ImV 
(lb-force)  a bias  of  12.  5 vdc  was  needed.  The  load  cell  bridge  had  an  offset 
of  17.  9 mV  with  no  load  applied. 

The  procedure  used  to  arrive  at  a value  of  p^^  began  by  taking  several  photo- 
graphs of  the  interferometer  output  for  a range  of  test  stand  pressure  levels. 
The  angle  of  the  tilting  of  the  vertical  fringes  was  measured  on  each  photo- 
graph. From  a measurement  of  horizontal  spacing  of  the  fringes  on  several 
photographs  an  average  fringe  spacing  d'  was  found.  Since  this  value  of 
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average  spacing  includes  a magnification  factor  M due  to  the  optical  system, 
the  true  value  of  d referred  to  the  sample  as  d = d'/M. 

A graph  was  made  of  vs  f and  a straight  line  fit  to  the  points  was  done 

d 

by  the  least  squares  approximation  in  the  linear  form 


tan  9 
d 


mf 


The  strain  optic  coefficient  p^^  is  proportional  to  the  slope,  m,  of  the  line 


ij 


T2^  "" 


(54) 


The  sign  of  p^^  was  determined  by  tilting  the  reference  mirror.  If  an  inward 
tilt  of  the  reference  mirror  (shortening  the  optical  path  length  of  the  top  of 
the  reference  beam  while  lengthening  the  optical  path  of  the  bottom  of  the 
reference  beam)  produced  the  same  tilt  on  the  fringe  lines  as  the  applied 
stress  had,  then  p^^^  was  positive,  otherwise  Pj^j  was  negative. 


Results  and  Discussion 


Three  different  materials  were  measured.  The  red  Corning  CS-2-60  filter 
glass  and  the  Lucite  PMS-1  samples  are  isotropic  and  only  p^j^  and  p^2 
values  are  required.  For  single  crystal  LiNbOg,  p^j,  P22*  P12  mea- 
sured using  an  X-cut  and  a Z-cut  wafer.  Other  quantities  required  wafers 
of  different  cut  and  were  not  available  at  the  time  of  this  measurement. 


Corning  CS-2-60  Glass  --  Typical  photographs  of  interference  fringe  pat- 
tern with  and  without  load,  as  taken  with  the  interferometer,  are  shown  in 
Figure  37(a)  and  (b)  respectively.  The  value  of  tan  ft/d  were  measured  for 
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I 
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(a) 


(b) 

Figure  37.  Photograph  of  fringe  lines:  (a)  with  no  load  and  (b)  under 
a load  of  26,  1 pounds,  measured  on  Corning  CS-2-60 
glass  sample. 
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various  forces  used  from  photographs  such  as  these  and  the  results  are  plot- 
ted in  Figure  38.  The  slope  of  the  least-squares  fitted  line  yields  the  photo- 
elastic coefficients.  The  final  results  are: 


p^j  = 0.  06  for  light  polarized  parallel  to  stress  direction 

and 

Pl2  = 0.  20  for  light  polarized  perpendicular  to  stress  direction 

The  Young's  modulus  used  for  the  above  calculation  are  assumed^ ^ to  be 
1.  02  X lo"^  psi. 

Lucite  PSM-1  --  Similar  measurements  were  made  on  Incite  and  the  results 
are  shown  in  Figure  39.  We  obtain 

Pjl  =0.  14 

p^2  “ 0*  25 

These  results  compare  well  with  the  manufacturer  specification  of  p^  ^ = 

0.  126  and  p^2  “ 0-  28.  It  is  interesting  to  note  that,  because  of  the  large 
photoelastic  effect,  the  sample  under  stress  causes  the  beam  to  deflect 
like  a prism  in  the  vertical  direction.  It  is  necessary  at  times  to  readjust 
the  reference  mirror  to  compensate  for  this  tilting  effect. 

LiNbOs  --  Because  of  the  limited  sample  size,  it  is  not  possible  to  com- 
pletely el-.ainate  the  fringing  stress  field,  therefore,  the  results  can  only 
be  used  as  a preliminary  reference  in  this  case.  As  mentioned  earlier,  two 
samples,  one  Z-cut  and  one  X-cut  were  used  for  this  measurement.  The 
data  are  shown  in  Figure  40.  The  results  as  well  as  the  published  data^^ 
are  given  in  the  following  tabulation: 
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Figure  38.  Photoelastic  coefficients  and  p^2  Corning  CS-2-60  glass 
as  measured  from  tan  0/d  vs  applied  force  plot. 
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Photoelastic 

Coefficient 

Our 

Results 

Published 

Results 

Pll 

-.  015 

-0.  02 

P22 

-.  015 

-0.  02 

P12 

0.  05 

0.  OH 

P31 

. 05 

0.  17 

In  conclusion,  we  have  performed  measurements  of  the  photoelastic  effect  on 
a number  of  samples  using  a Michelson  interferrometer.  The  technique  cjf 
measurement  involves  measuring  the  tilting  angle  of  the  reference  fringe 
lines  when  the  sample  is  under  load.  Thus  precise  knowledge  of  the  beam 
location,  and  mirror  tilting  are  not  necessary.  The  measurement  results  on 
samples  of  suitable  dimensions,  such  that  the  fringing  stress  field  due  to  the 
four-point  loading  arrangement  is  negligible,  are  reasonably  accurate  as  is 
demonstrated  by  the  good  comparison  between  the  published  data  and  our 
results  with  Lucite  PSM-1.  On  the  other  hand,  when  the  sample  dimension  is 
limited,  such  as  in  the  case  of  LiNbOg,  the  fringing  stress  field  introduced 
certain  error  in  the  measurements,  though  the  identical  results  obtained  for 
Pjj  and  Pj2  in  LiNb03  indicate  our  measurements  are  self-consistent.  For 
the  CS-2-60  red  filter  glass,  it  is  surprising  to  find  that  p^^^^  is  much  lower 
than  pj2*  ^nt  a p^2  value  of  0.  2 is  comparable  to  the  value  given  for  similar 
materials^. 


SKCTION  V 
CONCLUSIONS 


At  the  inception  of  the  research  investigation  on  the  strain-induced  wave- 
guiding  in  optical  material  for  integrated  optics  application,  we  have  outlined 
the  objectives  of  this  effort  as  follows; 

• To  determine  the  usefulness  of  the  strain-induced  effect  for 
integrated  optics  applications. 

• To  obtain  basic  understanding  and  perform  measurements  of  the 
effects  of  strain  on  the  various  physical  phenomena  in  optical 
materials  useful  for  integrated  optics  applications. 

• To  develop  device  concepts  using  selected  materials  and  methods. 
The  use  of  laser  beam  heating  for  strain-induced  integrated  optic 
circuits  will  be  emphasized. 

These  objectives  have  guided  our  research  effort  during  this  period  of  in- 
vestigation. Although  we  have  only  completed  phase  I and  the  initial  part  of 
phase  II  of  this  research  work,  substantial  progress  has  been  made  in  ful- 
filling these  objectives.  In  particular,  we  have  completed  all  tasks  originally 
outlined.  The  task  goals  were: 

• To  develop  techniques  of  inducing  controlled  residual  strain  for 
optical  waveguiding. 

• To  measure  the  effect  of  strain  on  selected  materials  and  to 
correlate  them  to  the  properties  relevant  to  optical  waveguiding. 
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• To  produce  strain -induced  planar  waveguides  or  channel  wave- 
guides and  obtain  the  waveguide  properties. 

• To  fabricate  simple  channel  waveguide  circuits  by  laser  heating. 

• To  construct  a microprocessor-controlled  facility  for  automated 
fabrication  of  laser-induced  channel  waveguides. 

In  this  report,  we  have  presented  the  research  results  and  progress  made 
toward  these  task  goals.  To  summarize,  we  have  achieved  the  following; 

• Demonstrated  laser  heating  strain-induced  (LHSI)  channel  wave- 
guides in  glass; 

~ Techniques  for  the  preparation  of  LHSI  channel  waveguide 
have  been  developed.  This  includes  the  determination  of 
the  laser  heating  temperature  profile,  the  relationship  of 
laser  power  and  scanning  speed  to  channel  waveguide 
dimensions  and  techniques  of  annealing  to  stabilize  the 
waveguide. 

LHSI  channel  waveguide  properties  have  been  measured  and 
analyzed.  This  includes;  the  waveguide  modes,  losses, 
input  coupling  angle  dispersion  and  refractive  index  profile. 

A model  for  understanding  how  the  LHSI  waveguides  are 
formed  has  been  developed.  The  strain  and  fictive  tempera- 
ture change  have  been  analyzed. 
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• Developed  fabrication  techniques  for  LHSI  channel  waveguide 
microoptic  circuits: 

Constructed  and  tested  a microcomputer  controlled  laser 
beam  exposure  system  (LBES)  for  automated  fabrication 
of  LHSI  waveguides. 

Fabricated  LHSI  waveguides  exceeding  2"  in  length. 

Fabricated  LHSI  waveguide  bends  and  gaps. 

Fabricated  LHSI  waveguide  directional  couplers. 

Constructed  and  tested  a fiber  optic  data  link  for  the  micro- 
processor. Controlled  two-dimensional  translating  stage 
for  LHSI  waveguide  fabrication. 

• Characterized  the  residual  or  induced  strain  in  materials  useful 
for  integrated  optics  applications: 

A polishing  procedure  to  introduce  a controlled  surface 
strain  in  LiNbO^  has  been  defined.  Both  grinding  and 
polishing  techniques  have  been  developed  to  produce 
residual  strain. 

X-ray  characterization  techniques  for  the  determination  of 
the  polishing-induced  strain  have  been  developed.  Measure- 
ment of  strain  using  Bond  techniques  has  been  obtained. 

- The  comparison  of  the  lattice  parameter  changes  in  LiNbO^ 
by  p>olishing-induced  and  Ti  indiffused  modifications  has  been 
completed. 
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A new  interferrometer  for  the  measurements  of  the  photo- 
elastic effects  in  material  medium  was  constructed. 

In  conclusion,  we  have  developed  the  laser  beam  heating  technique  for  the  fab- 
rication of  multimode  waveguide  circuits  automated  by  microprocessor  control. 
We  have  advanced  the  understanding  of  the  effect  of  strain  by  laser  heating 
and  optical  polishing.  The  major  goals  originally  set  forth  have  been  achieved. 
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